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^ ■ ABSTRACT 

(-H ' Context. Galactic planetary nebulae are observed with a wide variety of kinematic properties, spatial distribution, chemical compo- 

sition and morphologies, comprising members of the dominant stellar populations of our Galaxy. Due to their broad astrophysical 
I ■ interest, a proper characterization of these populations is of major importance. 

O ' Aims. In this paper we present a re-analysis of the criteria used to characterize the Peimbert classes I, Ila, lib. III and IV, through a 

statistical study of a large sample of planetary nebulae previously classified according to these groups. In the original classification, it 
C/3 , is usual to find planetary nebulae that cannot be associated with a single type; these most likely have dubious classifications into two 

or three types. Statistical methods can greatly contribute in providing a better characterization of planetary nebulae groups. 
Methods. We use the Bayes Theorem to calculate the posterior probabilities for an object to be member of each of the types I, Ila, lib, 
III and IV. This calculation is particularly important for planetary nebulae that are ambiguously classified in the traditional method. 
The posterior probabilities are defined from the probability density function of classificatory parameters of a well-defined sample, 
composed only by planetary nebulae unambiguously fitted into the Peimbert types. Because the probabilities depend on the available 
observational data, they are conditional probabilities, and, as new observational data are added to the sample, the classification of the 
nebula can be improved, to take into account this new information. 

Results. This method differs from the original classificatory scheme, because it provides a quantitative result of the representativity 
of the object within its group. Also, through the use of marginal distributions it is possible to extend the Peimbert classification even 
^\ • to those objects for which only a few classificatory parameters are known. 

' Conclusions. We found that ambiguities in the classification of planetary nebulae into the Peimbert types, should be associated to 

, difliculties in defining sharp boundaries for the progenitor star mass for each of these types. Those can be at least partially explained 

. by real overlaps of some of the parameters that characterize the different stellar populations. Those results suggest the need of a larger 

number of classificatory parameters for a reliable physical classification of planetary nebulae. 
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1. Introduction served with a wide variety of spatial and kinematic properties, 

chemical composition, and morphology, comprising a mixed 

Planetary nebulae (PNe) result from the ejection of the external g^^up of objects ranging from Population I (slow moving disk 

envelope of low and intermediate mass stars during their evolu- objects) to extreme Population II (fast moving objects from 

tion from the asymptotic giant branch (AGB) to the white dwarf j^e bulge and the halo). On account of this, PNe are objects 

stage. Roughly speaking, the nebula will be created by the in- j^^gg astrophysical interest. Their application is multidisci- 

teraction of the slow AGB wind with the new fast stellar wind pjin^ry: they probe the evolution of low and intermediate mass 

from a post- AGB phase, which will sweep up the circumstellar ^^^.^^ ^^e Galactic nucleosynthesis history. Elements that are 

envelopes into a shell-like structure (Kwok et al. 1978; Kwok modified by the nucleosynthesis of the parent star allow us to 

1982, 1983; Kahn 1983 and Yolk & Kwok 1985). The energy investigate the composition of the interstellar medium at the time 

source of the nebula is the ultraviolet radiation from the hot cen- ^^en the star formed. PNe are useful tracers of Galactic kine- 

tral star, which gradually ionizes the shell, at the same time the j^^tics (Schneider & Terzian 1983; Durand et al. 1998; Maciel 

shell grows in mass as more AGB wind material is swept up. & ^^g^ 2005), Galactic chemical evolution (cf. Costa & Maciel 

Because PNe come from stars in a large range of main se- 2006), and are even useful tools in the investigation of some 

quence masses (0.8-8.0 Mq, cf. Peimbert 1990), they are ob- cosmological problems (see Balser et al. 1997, 1999, 2006, and 

references therein about the measurement of ^He abundance in 

Send offprint requests to: C. Quireza PNe and astrophysical implications). This richness of informa- 

* Table [T] and Table |6] are only available in electronic at the tion generally requires a proper selection of the objects to be 

CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.125.5) or via used in a specific analysis. By selecting PNe that closely reflect 
[http://cdsweb.u-strasbg.fr/Abstract.htmli 
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the properties of the interstellar medium out of which their cen- 
tral stars have been formed, we can study, for instance, the pres- 
ence of abundance gradients in the disk of the Galaxy. There is 
not an ideal method of performing such selection, since how rig- 
orous one is depends on the kind of study to be performed and of 
the goal of the analysis. However, it is important to understand 
the different groups (or populations) of PNe and, at least, some 
of their main characteristics. 

Among the different planetary nebula classification schemes 
found in the Hterature (Greig 1971, 1972; Kaler 1983; Heap & 
Augensen 1987), one of the most efficient is probably that pro- 
posed by Peimbert (1978), which suggested the classification of 
PNe into four types, namely I, II, III and IV (see also Peimbert & 
Serrano 1980; Peimbert & Torres-Peimbert 1983). Each of these 
groups should roughly correspond to a stellar population of the 
Galaxy. These groups also represent a sequence of mass inter- 
vals for the planetary nebula progenitor star, type I representing 
the youngest population, with the most massive progenitors, and 
type IV the oldest population, with the less massive progenitor 
stars. Later a type V was also defined to include bulge objects 
(Maciel 1989), which should constitute a distinct stellar popu- 
lation, but the selection criteria are different than those used to 
select the four original types. In order to characterize the types 
I to IV, Peimbert took into account, chemical, kinematic al and 
spatial properties of PNe. Such criteria were supplemented by 
Peimbert & Torres-Peimbert (1983), Faundez-Abans & Maciel 
(1987a) and Maciel (1989), and have been continuously used 
in Galactic evolution related issues, mainly in the study of the 
abundance gradients (Maciel & Faiindez-Abans 1985; Faiindez- 
Abans and Maciel 1987b; Maciel & Quireza 1999; Maciel 2002; 
Maciel et al. 2003, 2006; Perinotto & Morbidelh 2006). The 
classification, however, has its flaws: some criteria are poorly de- 
fined, and there is a tendency to find objects with different char- 
acteristics in the same class; moreover, several objects simply do 
not fit in any of the existing classes, which suggests the need for 
a larger number of classificatory parameters, and/or more rigor- 
ous classification criteria. 

In this paper we present a re-analysis of the criteria used to 
characterize the Peimbert classes, through a statistical study of 
a large sample of PNe previously classified into the types I to 
IV. The process of selection of the fundamental parameters is 
described in ^ In |3]we describe the classification criteria by 
Peimbert (1978) and we present the classification of our PNe ac- 
cording to such criteria. In ^we present our statistical analysis, 
which provides us information about how well each of the ob- 
jects accommodates in a given class. We discuss the results and 
implications of our analysis in ^ In ^we summarize our main 
results. 



the Strasbourg-ESO Catalogue of Galactic Planetary Nebulae 
(Acker et al. 1992), plus new observed PNe, not present in the 
Acker et al. catalogue (Kohoutek 1994; Beaulieu et al. 1999; 
Escudero & Costa 2001). The resulting sample has a total of 
476 objects. The fundamental parameters we looked for were 
the helium, nitrogen and oxygen abundances relative to hydro- 
gen, the heliocentric distance, d^^i, and the radial velocity rela- 
tive to the Local Standard of Rest, Vlsr- Additionally, we have 
compiled the flux at 5 GHz, S 5 ghz, and the angular diameter of 
the PNe, ©diam, which have been used to exclude bulge objects, 
as will be discussed in § |3] The adopted parameters are listed in 
Table [T] which is available only in electronic form. The process 
of selection of these parameters is described in the subsections 
ahead. In the table, for each planetary nebula, we also list the 
common name, the PNG number (Acker et al. 1992), and cal- 
culated properties, namely the Galactocentric distance, Rcai, the 
height in relation to the Galactic plane, z, and the peculiar ve- 
locity. Ay. Explicit formulae for these quantities are given later 
on in ^3.1! The resulting Peimbert classes are listed in the last 
column. Those quantities will be commented in § [3] For units, 
the distances are given in kpc, the angular size in arcsec, and the 
flux in mJy. The oxygen, and nitrogen abundances are given as 
log(X/H) + 12, as usual. The corresponding references for each 
of these quantities are provided in the table. Most of the Galactic 
coordinate values come from the Acker et al. (1992) catalogue. 

2.1. Elemental abundance 

For each planetary nebula, we have tried to compile all nebular 
abundance measurements available in the literature. The sam- 
ple should be as complete as we could make it up to 2004. 
Additional abundances published after 2004 have been occasion- 
ally included, such as the ISO (Infrared Space Observatory) re- 
sults from Pottasch & Bernard-Salas (2006). For each object, 
we have kept only those abundance determinations that agree 
within the adopted uncertainties, which we estimate to be in the 
range 0. 1-0.4 dex for most objects. The adopted abundances are 
the average of all the remaining abundances for a given PNe. 
Abundances from Costa et al. (1996, 1997, 2004), Cuisinier et 
al. (2000), Escudero & Costa (2001), and Escudero et al. (2004) 
were preferred to be used individually or averaged. We did so 
because such studies used the same instruments, similar obser- 
vational techniques and data analysis procedures. Pottasch & 
Bernard-Salas (2006) were also included in this last group, since 
their abundance measurements should be one of the most accu- 
rate available in the literature up to now. 



2. Planetary nebulae data base 

Ideally, when analyzing a number of objects, it is desirable to use 
a homogeneous sample, with data obtained using similar tech- 
niques and self-consistent measurements, which results in very 
accurate information. Unfortunately, samples with precise and 
self-consistent values of abundance, distances, and other funda- 
mental parameters of PNe are still not large enough for statistical 
studies. We have concentrated efforts in generating a large sam- 
ple of PNe of different properties, so we can characterize the var- 
ious populations/groups of PNe, belonging to the different struc- 
tures of our Galaxy. In order to do that, we have compiled the 
nebular properties from selected sources from the literature. 

Our sample consists of PNe for which abundances have been 
determined in the literature. Most of it consists of objects from 
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+ e(X) = log(X/H) + 12. 

* Radial velocities from the catalog of Durand et al. (1998). For every object, the heliocentric radial velocities have been converted to the Local Standard of Rest radial velocities (see text). 
References of abundances. (3) Costa et al. 2004; (21) Perinotto 1991; (22) Perinotto et al. 1994; . . . 

References of distances. (4) Costa et al. 2004 (see references in); (7) Harris et al. 1997; (10) Maciel 1984; (17) Pottasch 1996; . . . 
References of the angular diameters. (7) Cahn & Kaler 1971; (18) Perek & Kohoutek 1967; (27) Zhang 1995 (see references in); . . . 
References of the 5 GHz flux densities. (9) Zhang 1995 (see references in); . . . 
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2.2. Heliocentric distances 

The range in distances as determined by different authors can be 
very large, and choosing a good distance estimate is a difficult 
task. Only a small number of PNe have accurate individual esti- 
mates of their distances. For the vast majority of PNe, the only 
distances available are those obtained from statistical methods. 
Here, whenever possible, we use individual distances over sta- 
tistical estimates. Those include trigonometric parallax (Harris 
et al. 1997), Hipparcos parallaxes (Pottasch 1997; Pottasch & 
Acker 1998), spectroscopic parallax, radial expansion distances 
measured by VLA or by Hubble Space Telescope (Hajian et al. 
1995; Terzian 1997) and local extinction (Martin 1994). We have 
chosen as far as possible to use distances estimated from paral- 
lax or expansion method, because in both cases the distances are 
independent of any (assumed) property of the nebula or its cen- 
tral star. In Table [1] individual distances represent about 17% 
of the sample. Among statistical distances the catalog of Maciel 
(1984) represent almost half of the distance sample, and is com- 
plemented by Cahn et al. (1992) and Zhang (1995). Cahn et al.'s 
distances also substituted lower and upper limits from Maciel's 
catalog. For the bulge and anticenter PNe compiled respectively 
from Escudero et al. (2004) and Costa et al. (2004), original ref- 
erences for the distances are Schneider & Buckley (1996), van 
de Steene & Zijlsti-a (1995), and Amnuel et al. (1984). We could 
not find distances for 71 objects in our sample. 

2.3. Radial velocity relative to the Local Standard of Rest 

We have used heliocentric radial velocities from the catalog of 
Durand et al. (1998), which results from the selection of all 
known Galactic PNe radial velocities up to 1998. For a given 
object in their catalog, radial velocities are weighted averages of 
each existing velocity by the inverse of the square of its asso- 
ciated error. According to the authors, 90% of the total sample 
have accuracies better than 20 km s"^ 

The heliocentric radial velocities, Vhd, have been converted 
to the Local Standard of Rest radial velocities using the basic so- 
lar motion (u, v, w) - (-9, 11,6) kms"', as given by Mihalas & 
Binney (1981). For aright hand coordinate system the coiTection 
is 

Vlsr = Vhei + (9 COS COS Z + 1 1 cosfe sinZ + 6 sinb) kms"' , (1) 

were I and b are the Galactic coordinates of the objects. We have 
radial velocities for 428 PNe in our sample. 

2.4. Angular diameter 

Whenever possible, we have used radio continuum measure- 
ments over optical ones. We also have tried to give priority 
to measurements made by a same group of researchers, which 
generally use the same observational techniques (for instance, 
Aaquist & Kwok 1990, 1991; Kwok & Aaquist 1993). Those 
choices have been taken as an attempt to make the sample less 
heterogeneous. In Table [T] we list angular sizes for 448 objects. 
For elongated nebulae, a geometric mean of the semi-major and 
semi-minor axes was used whenever necessary. 

2.5. Radio continuum flux density at 5 GHz 

Most of the radio fluxes have been adopted from Zhang (1995), 
which used data obtained with the Very Large Array (VLA) 
(Aaquist & Kwok 1990; Zijlstra et al. 1989; Gathier et al. 
1983), plus some single-dish measurements (Milne & Aller 



1975; Milne 1979) and measurements tabulated by Pottasch 
(1984) (see Zhang & Kwok 1993). We list fluxes for 388 ob- 
jects. 

Angular sizes and flux densities have been used with the sole 
purpose of selecting bulge objects in our sample (§ O. These 
quantities have not been used in our statistical analysis. 

3. Pre-classification of planetary nebulae 

3.1. Peimbert types and classification 

Peimbert (1978) has divided the Galactic PNe into four types 
on the basis of their chemical, spatial and kinematic properties. 
Ideally, each of these groups should correspond to distinct inter- 
vals of mass of the progenitor star, whose evolution should af- 
fect differently the chemical composition of the ejected planetary 
nebula envelope (Calvet & Peimbert 1983; Peimbert & Serrano 
1980; Peimbert & Torres-Peimbert 1983). 

Type I PNe are He and N rich, presenting extremely filamen- 
tary structure. They have been defined by Peimbert (1978) as 
PNe with He/H > 0.14 or log (N/O) > 0, but later on Peimbert 
& Torres-Peimbert (1983) have relaxed this definition by includ- 
ing those objects with He/H > 0.125 or log (N/O) > -0.30. 
This excess in He and N can be eventually measured as a N/O 
excess in the nebula (Peimbert 1990), and has two possible ex- 
planations: i) N contamination in the outer atmosphere of the 
planetary nebula progenitor star, caused by dredge up processes. 
Such N would be produced at the expenses of C and O during 
H burning by the CNO cycle; and ii) The progenitor star would 
have been formed more recently from a medium richer in He 
and heavy elements. Both causes suggest that type I PNe result 
from the most massive progenitor stars in the intermediary mass 
range, being the youngest of all disk PNe. They are associated 
with the Galactic thin disk (closer to the Galactic plane), whose 
scale height is, here, assumed to be about 300 pc (Gilmore & 
Reid 1983, Reddy et al. 2006). Because thin disk stars orbit 
the Galactic center in nearly circular orbits, type I PNe should 
present low velocity dispersion (Dutra & Maciel 1990; Maciel 
& Dutra 1992). 

Type II PNe belong to the intermediate disk population. Still 
associated with the thin disk, they have approximately circular 
orbits. Members of this group are generally older than type I 
PNe, having been formed at a time when the interstellar medium 
was more metal-poor in heavy elements. Therefore, they may 
present slight underabundances of elements such as O, S and Ne. 
Moreover, their low mass prevents strong He or N enrichments. 
According to Peimbert (1978), most existing PNe are probably 
of type II, having He/H < 0.125 and log(N/0) < -0.3. Later on, 
Faundez-Abans & Maciel (1987a), further subdivided the type II 
PNe into types Ila and lib, according to their N abundance (types 
lib having no N enrichment and types Ila having some N enrich- 
ment intermediate between type I and type lib). In such scheme, 
type Ila have been formed from progenitor stars near the high 
mass bracket. Because of that, type Ila PNe may present slightly 
enriched abundances of O, S, Ne, and Ar, and some definite en- 
richment in He and N. 

Here, as in Maciel & Quireza (1999), we consider a slightly 
more rigid criteria to sample PNe in the thin disk: we con- 
sider as type I only those objects for which both conditions 
He/H > 0.125 and log(N/0) > -0.30 ai-e satisfied. As a conse- 
quence, objects for which only one of these conditions is satis- 
fied have been considered as type Ila (Faundez-Abans & Maciel 
1987a), so that they are in fact type II objects. This procedure 
should assure a correct sampling of the nebulae associated with 
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the more massive progenitor stars, in the high end of the mass in- 
terval that characterizes intermediate mass stars. In this interval, 
which corresponds to type I planetary nebula, the nebular oxy- 
gen abundance may be affected by ON cycling in the progenitor 
stars (Henry 1990; Maciel 1992; Perinotto & Morbidelli 2006). 

Type III PNe have been defined by Peimbert (1978) as those 
objects with |Ay| > 60 km s"' which do not belong to the halo 
population. They are high velocity objects from the thick disk 
population. Such PNe were ejected from stars, which are gen- 
erally older than most thin disk stars. As most of the thick disk 
stars they are on moderately elliptical orbits that typically reach 
higher distances from the Galactic plane. Type III nebulae should 
also present underabundances of heavy elements, compared to 
thin disk stars. Here we assume that the thick disk have a scale 
height of 1450 pc (Gilmore & Reid 1983, Reddy et al. 2006). 

Finally, type IV PNe are halo objects. Only a few observed 
objects belong to this class, identified by the low heavy element 
abundances and the large deviations from disk kinematics (high 
velocity with respect to the Sun). These nebulae are remnants of 
very low mass halo stars (a rather uniform old and metal-poor 
population). They do not show excess in helium, but they appear 
to have a slight deficiency of this element relative to the other 
PNe. 

There is still a fifth group named type V, which consists 
of the nebulae found in the Galactic center (Maciel 1989). 
According to Stasihska et al. (1992), this group is composed 
by 700 known PNe. Bulge objects have a distinct evolutionary 
history than the disk and halo and a variety of chemical compo- 
sition and stellar masses might be related to this component of 
the Galaxy. For this reason, type V nebulae cannot be defined 
in terms of the same criteria used to define the types I to IV. 
Ignoring bulge PNe as a distinct population, however, leads to 
larger uncertainty in the classification of the other types, as will 
be discussed separately in § 13.21 Here we used a method sim- 
ilar to that used by Pottasch & Acker (1989), Stasinska et al. 
(1991) and Zhang (1995) to select bulge objects: i) angular ra- 
dius smaller than 10 arcsec; ii) Galactic longitude (I) and latitude 
(b) within 10° of the Galactic center (which places the nebulae 
1.32 kpc above or below the Galactic plane, if Rq - 7.6 kpc, see 
Maciel 1993), and iii) flux at 5 GHz less than 100 mJy (Pottasch 
1990). Because the criteria used to select type V are not common 
to those adopted for the types I to IV, this group cannot be used 
in our statistical analysis (§|4]i. 

For types I to IV, the limits between each group are defined 
by four quantities: the He/H and log(N/0) abundance ratios, the 
absolute values of the height above the Galactic disk, |z|, and of 
the peculiar velocity, |Ay|, of each object, following that order 
of importance. 

The height above the Galactic disk, z is given by 



where, /?Gai is the distance to the Galactic center projected in the 
Galactic plane. 



z - dhei sin b , 



(2) 



where t/hei is the heliocentric distance of the object. The pecu- 
liar velocity is basically the difference between the measured ra- 
dial velocity relative to the Local Standard of Rest, Vlsr, and 
the radial velocity obtained from the Galactic rotation curve (as- 
suming the object has a circular orbit). The planetary nebula ra- 
dial peculiar velocity is calculated from the equation (Maciel & 
Dutra 1992) 



[9 9 
^0 + (dhei COS b) -2Rq <ihei cos b cos ZJ 



(4) 



0(^Gai) is the rotation velocity at the nebula galactocentric po- 
sition calculated for a rotation curve. We adopted the rotation 
curve given by Maciel & Dutra (1992) for type I + II PNe. 



0(^Gai) = 314.1356- 20.5234 Rg^i + 0.9855 rI^^ ., 



(5) 



^Gal ^0 



sin Zcosi» , 



(3) 



For the rotation velocity in the Sun position we adopted 0© = 
185 kms"' (Rohlfsetal. 1986). 

The criteria for classification of the four types proposed by 
Peimbert (1978), after some reassessments (Peimbert & Torres- 
Peimbert 1983; Faiindez-Abans & Maciel 1987a; Maciel 1989; 
Maciel & Koppen 1994), are shown in Table |2] The last three 
columns give, respectively, mean estimates for the progenitor 
mass and the age intervals corresponding to each type, and the 
expected Galactic population to which the planetary nebula be- 
longs. For completeness, the type V PNe were also listed in the 
table. 

It should be noted that the masses and ages given in Table |2] 
are based on the original Peimbert criteria, and do not take into 
account some more recent attempts to estimate these parameters 
on the basis of individual methods, as in Maciel et al. (2003, 
2005, 2006). Mass intervals in Table|2]should be used as a guide 
only, since large uncertainties are associated with estimates of 
progenitor and central star masses. Besides, recent results from 
Marigo (2007) suggest that hot bottom burning (supposed to 
be responsible for an overabundance of He and N in the stel- 
lar outer layers) may be weakened or even prevented by the 
third dredge-up process during the early stages of the thermally- 
pulsating AGB phase in stars within « 3.0-4.0 Mq and metal- 
licities Z 0.001, if effects of variable molecular opacities are 
taken into account. This could reduce the mass range of type I 
PNe to masses higher than x 3.0-4.0 Mq. 

The resulting classification, following the criteria displayed 
in Table|2]is shown in the last column of the Table[T] In this table, 
34% of the whole sample are type II objects (24% Ila and 10% 
lib), 28% are type V and 9% are type I. Only 3 objects were clas- 
sified as type III and 7 objects were classified as type IV (each 
group representing 1 % of the whole sample). Approximately 4% 
of the sample (without He/H abundance) were not classified. The 
rest of the sample (approximately 23%) consists of objects that 
could not be fitted in one class, most being between types II, III 
and IV. 

In spite of having those criteria, the classification process is 
not straightforward since most nebulae do not satisfy all the nec- 
essary criteria to fit undoubtedly into a given class. Most likely, 
they can be classified into two or three types or not even being 
classified at all. One example is NGC 6807. For this planetary 
nebula He/H = 0.110, log(N/0) = -0.60, |z| = 0.66 kpc, and 
|Ay| - 162 kms"'. Its position indicates that this object can 
be a type lib nebula. The high peculiar velocity, however, indi- 
cates that this object may be a type III or IV. Because there is 
no limit in abundance among these three groups, the planetary 
nebula cannot be classified in one single class. 

In an attempt to evaluate how accurately each object fits in 
its group, some flags were added to the resulting classification. 
We use an asterisk "*" to indicate absence of log(N/0) mea- 
surement (in this case the object cannot have an accurate classi- 
fication). Letter "A" indicates that one of the conditions between 
angular radius smaller than 10 arcsec and flux at 5 GHz less than 
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Table 2. Classification criteria for Peimbert classes 



Type 




He/H" 


log(N/0)" 


log(N/H) + 12 


(knc) 


|AV| 
(kms"') 




Age 
(Gv) 


Location 


I 




> 0.125 


> -0.30 




« 1 (< 0.3) 


< 60 


2.4-8. 0'^ 


0-2'' 


thin disk 


II 


Ila 
Ila 
lib 


> 0.125 

< 0.125 

< 0.125 


< -0.30 
> -0.60 

< -0.60 


A IV IV 

OO 00 00 

bob 


< 1 


< 60 


1.2-2.4"^^ 


4-6'' 


thin disk 


III 










> 1(< 1.45) 


> 60 


1.0-1.2'^ 


8-10'' 


thick disk 


IV 










» 1 (> 1.45) 


> 100 


o.s-i.o-^ 


> 10 


halo 


V 










< 1.3 


large range 


large range 


large range 


bulge 



^ No abundance criteria were originally defined to the types III, IV and V. Average abundances can be calculated using data from Chiappini & 
Maciel (1994) for types III and IV, and data from Escudero & Costa (2001), Cuisinier et al. (2000) and Escudero et al. (2004) for type V. Mean 
values of He/H and log(N/0) are respectively 0.099 and -0.68 for type III, 0. 104 and -0.67 for type IV and 0. 1 10 and -0.60 for type V PNe. 

^ According to the original criteria, type I nebulae have |z| <c 1 and type III |z| > 1 . Here, quantitative limits for |z| were defined assuming the 
scale height of 0.3 kpc to the thin disk and 1.45 kpc to the thick disk (Gilmore & Reid 1983). Criteria for type V are defined using = 7.6 kpc 
and assuming that the object should be within 10° from the Galactic center. 
Maciel (1992) 

'' Maciel &K6ppen (1994) 



100 mJy is not obeyed (applies only for type V nebula). Letter 
"B" indicates that the condition for |z| is not obeyed. Letter "C" 
indicates that the condition for |Ay| is not obeyed. Letter "D" 
indicates absence of |z| and/or |A V| (and/or angular radius and/or 
flux for the type V nebulae). Letter "E" indicates that more than 
one among the conditions "A", "B", "C" and "D" are applied 
(for types I, Ila and V). Finally, a "F" indicates a perfect clas- 
sification: the source has all characteristics of the type. In the 
example cited in the last paragraph, the resulting classification 
of NGC 6807 would be IIbC/(IIIAV)B. 

3.2. Reasons for a posterior classification 

The merit of the Peimbert classification is in the attempt of sam- 
pling PNe in different populations (halo, thin disk, thick disk). 
Naively, the sequence of types I-II-III-IV could be translated 
into a sequence of decreasing mass intervals for the planetary 
nebula progenitor stars in the main sequence, in the sense that 
type I have high mass progenitors, type II intermediate mass, 
and type III and IV low and very low masses, respectively. In 
reality, although distinct mass intervals are usually provided for 
each of these types, these are only mean estimates, and well de- 
fined boundaries between them cannot be established from the 
four criteria, He/H, N/O, |z| and |Ay| only. Peimbert himself, in 
Peimbert & Carigi (1998), affirmed that it is not possible to have 
a sharp mass boundary between types II and III based on dy- 
namical arguments only or a sharp mass boundary based on the 
observed abundances, either. Three causes are provided: i) ob- 
servational errors in the abundances determinations, ii) displace- 
ments of the older PNe from the original Galactocentric distance 
at which the progenitor star originated, and iii) a real scatter in 
the chemical abundances of the interstellar medium at a given 
Galactocentric distance. 

Besides those reasons, additional factors may contribute to 
difficulties in the classification of PNe into types I to IV: 

1 . Properties that are common to distinct stellar populations: 
although the mean properties of the dominant stellar popu- 
lations of the different components of the Galaxy are rea- 
sonably well defined, some characteristics may be common 



to them (Wyse & Gilmore 2005). This problem is reflected 
in the He/H and N/O abundance ratios for types lib. III 
and IV PNe. All these types are low metallicity objects, 
so they are likely to have low He and O abundances. At 
first a low N abundance would also be expected both be- 
cause: i) they come from an old stellar population, proba- 
bly formed from a nitrogen poor interstellar medium; and 
ii) they are remnants of low mass stars so the nebula should 
not be self-contaminated by He and N from the stellar nu- 
cleosynthesis. Based on this kind of reasoning, we should 
find progressively lower He, O and N abundances as we go 
from type lib to type IV PNe. The interpretation in terms 
of the N/O ratio, however, is not as simple. From the limit 
of log(N/0) for type lib in table |2] one may expect to find 
log(N/0) < -0.60 for both thick and halo PNe, but obser- 
vations show that, at least in regard to the halo population, 
this is not true. Howard et al. (1997) measured the chemi- 
cal composition of nine halo PNe, finding that most objects 
show enhanced log(N/0) abundance raticQ. The spread in 
log(N/0) is much larger than can be accounted by uncer- 
tainties alone. Unfortunately, the small number of type IV 
PNe with reliable abundances and distances somewhat hin- 
ders their application in statistical studies. From compiled 
data up to the end of 2001, Stasinska (2004) mentions a total 
of only 20 known halo PNe (Wyse & Gilmore 2005; Costa 
& Maciel 2006). The large scatter in the N/O abundance ra- 
tio in metal-poor stars is also suggested by high precision 
abundance data for Galactic halo stars (Spite et al. 2005; 
Israelian et al. 2004). In these papers, the scatter in the N/O 
abundance ratio is also much larger than their quoted error 
bars. A model able to predict high N/O abundance ratios at 
low metallicities was presented in Chiappini et al. (2005), 
where the explanation for the high nitrogen abundance was 
based on an increase of the rotational velocity in very metal- 
poor stars (Meynet & Maeder 2002). More recent models 
by Chiappini et al. (2006a, b), considering N yields for stars 
above 20 Mq at metalUcity Z = 10"^ (Hirschi 2007), sug- 



' Values range from -0.91 to -1-0.17; only 3 out of the 9 PNe have 
log(N/0) < -0.60. All nine objects exhibit subsolar 0/H. 
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gest that at log(0/H) + 12 < 7 (or [Fe/H] < -3), the main 
contribution for the larger amounts of N would come from 
fast rotating massive stars. After log(0/H) + 12 ~ 7, where 
AGB contribution would start to be effective, intermediate 
mass stars of high rotational velocities could also produce 
large amounts of nitrogen. In this framework, the scatter in 
N/O abundance ratios would reflect the distribution of stel- 
lar rotational velocities as function of metallicity. In spite of 
the uncertainties associated to the observed data, this picture 
suggests that the interstellar medium from which the halo 
PNe progenitors have been formed had a significant scat- 
ter in N/O. Thus, even considering that the contribution of 
the halo progenitor star for the N enrichment should not be 
significant (once these stars are of very low mass) the use 
of PNe N/O abundance ratio as an accurate tracer of stellar 
populations is unclear. Besides, a real scatter in abundance, 
another complication arises if the position of the PNe was 
affected by the motion of the progenitor star away from its 
original galactocentric birth radius. This possibility is par- 
ticularly suitable for the older and high velocity PNe. As a 
consequence, it is possible that among PNe considered to be 
in the halo, some actually probably belong to an old disk 
population; also, samples of PNe considered as members of 
the Galactic disk may actually contain PNe of the halo pop- 
ulation that are found physicaUy in the same region as the 
disk. We may also find PNe of the disk or halo population 
that are found physicaUy in the same region as the Galactic 
bulge. 

2. Uncertainties in the classification parameters and the ab- 
sence of a homogeneous data sample: we verified that typi- 
cal differences in abundance measurements resulted mostly 
in changes in the classification between types Ila to I and 
lib to Ila. Nevertheless, the number of objects for which this 
happens is really small. A major source of uncertainty in the 
classification comes from the adopted distance. The use of 
a different distance scale does not have a significant effect 
over the classification of type II and I nebulae, for which the 
classification criteria are more dependent on chemical abun- 
dances. However, some uncertainties may arise in the classi- 
fication of types lib. III and IV, for which classification has 
a larger dependence on distance-related parameters. 

3. Objects located in the direction of the Galactic center: the 
classification of PNe seen towards lines of sight that are in 
the direction of the Galactic center has a larger degree of 
uncertainty than usual. This is due to the broad interval of 
possible values for the observed radial velocity. Because a 
variety of radial velocities are allowed to be observed in the 
direction of the center of our Galaxy, no limit can be estab- 
lished over the peculiar velocity, AV, and it cannot be used 
as a constraint anymore. As a consequence, instead of four 
parameters (namely He/H, log(N/0), |z|, |Ay|), only three 
criteria can be used to classify such objects, increasing un- 
certainties in the classification. This problem is significantly 
minimized when we consider PNe from the Galactic center 
as a distinct population (type V or bulge PNe), and exclude 
this groups from the sample previous to the classification of 
the PNe sample into the types I, Ila, lib, IIII and IV. 

4. Posterior classification of planetary nebulae 

When one of the classificatory parameters is not available to a 
given PN, it is usually not possible to associate this object to 

a single Peimbert type. A more rigorous procedure, based on 
the observational data available to the PNe, uses the posterior 
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probabilities for an object to be a member of a given type that 
is defined from the distribution of parameters of a well-defined 
sample. 

We consider an ideal sample for which all PNe are satisfac- 
torily distributed into the Peimbert types, each nebulae having 
classificatory parameters that satisfy all the conditions imposed 
by the classificatory criteria. The distribution of classificatory 
parameters (He/H, log N/O, |z| and |Ay|) of the population of ob- 
jects in each group may be estimated by a multivariate parameter 
distribution of a sample of such group. The probability of each 
planetary nebula in particular to be member of a given group 
is calculated by using the multivariate distribution of the group. 
This method differs from the original classificatory scheme, be- 
cause it provides a quantitative result of the representativity of 
the object within its group. Also, through the use of marginal 
distributions it is possible to extend the Peimbert classification 
even to those objects that do not have aU the necessary classifi- 
catory parameters. 

In addition to type V PNe, we have also excluded from the 
analysis those PNe having He/H < 0.060 and those having si- 
multaneously very high He/H (generally He/H > 0.180) and 
very low N/O abundance ratios. In the first case, measurements 
of He/H < 0.060 are Ukely to be errors since these values are 
lower than the primordial helium abundance (assuming 0.2384 
for the pre-Galactic helium abundance by mass, from Peimbert 
et al. 2002, which should correspond to an abundance by number 
of about 0.060). PNe that are observed to have very low He/H 
abundances should be of low-excitation, having some amount 
of neutral helium present in the nebula. Because neutral helium 
cannot be measured, the total helium abundance is obtained from 
the sum of the He^/H^ and He^^/H^ ionic abundance only. In 
those cases, the measured He/H abundance cannot be derived ac- 
curately, being underestimated in most cases. For the high He/H 
abundance ratio, we expect to see a trend of increasing N/O as 
He/H increases. This trend is commonly explained by dredge- 
up episodes, which occur in the red giant branch (RGB) and in 
the AGB phases: in particular, for the most massive stars ex- 
periencing hot-bottom burning during the thermal pulses on the 
AGB phase, the abundance of '"^N should be significantly af- 
fected, due to the conversion of C into N in the envelope. The 
hot-bottom burning yields higher N/O and lower C/O ratios. 
(Marigo 2001; Perinotto et al. 2004). Finally, it should be noted 
fliat, although He2-152 (He/H = 0.269 , log(N/0) = +0.60) 
and Ml-11 (He/H = 0.094 , log(N/0) = +0.77) do not fit in 
the group of PNe with simultaneously very high He/H and very 
low N/O abundance ratios, those PNe are outUers in the trend of 
increasing N/O versus He/H and for this reason they were also 
excluded. 

Let us consider x a four-dimensional vector composed of 
the four classificatory parameters of the K Peimbert types: x = 
(He/H, log N/O, |z|,|AV|). The probabiUty of finding a nebula 
having the set of observables x is: 



K 

p{x\D = YjP(.k\T)p{x\k,D, (6) 

where p{x\k, I) is the likelihood that one can observe x given 
that the nebula belongs to the group k and p(k\I) is the prior 
probability thai a given nebula belongs to the group k, when no 
other information is provided. The symbol I emphasizes that 
these probabilities depend on the state of the available back- 
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ground information. The prior probability obeys the normaliza- 
tion condition: 



(7) 



and p(k\T) > 0, Wk. We assume that p(x\k, T) belong to the same 
parametric family of distributions for all k, that is, the likeli- 
hood functions have the same profile and are distinguished in 
each group by a small number of parameters. We assume that 
pix\k, I) is a four-dimensional Gaussian, represented by 



pix\k,D^{27Tr^\^,\-''^ 



exp 



■^{x-fikfi:,'ix-fik) \ , (8) 



*-^He/H'^logN/0'^l 



fJ-^^Y^) is the vector of means, 



where fii^ 

is the covariance matrix corresponding to group k, and is its 
determinant. 

The marginal likelihood of a variable Xq may be found by 
integrating the p{x\k, I) over all other variables: 



p(x,\D = J] ■■■ pix\k,I)Y\dxi 

K 

= Y,p^k\i)p(x,\k,i), 



(9) 



k=\ 



where i stands for the four classificatory parameters used in the 
Peimbert classification. From Eq.|9]we see that the marginal like- 
lihood of Xq, of a population mixture, is composed by the mixing 
of the marginal likelihoods of Xq in each of the K groups consid- 
ered. Also, the marginal likelihood of any observable set y com- 
posed by a sub-group of variables of x is given by the integration 
of p{x\k, I) over the other variables which do not belong to y. 

Our goal is to estimate the probability that a nebula belongs 
to a given Peimbert type after considering all its available ob- 
servational data. According to Bayes Theorem, we can calculate 
these posterior probabilities as 



p{k\y,I) 



p{k\I)p(y\k,I) 

Z%iP(j\i)p(y\j,i) 



(10) 



where y is the set of known classificatory parameters for the neb- 
ula and p(y\k,T) is the marginal likelihood of observing the set 
_y in a nebula belonging to group k. Note that the posterior prob- 
abilities satisfy the equation ^^.^j p(k\y,T) - 1, Vj. As new ob- 
servational data are added to the sample, the classification of the 
nebula can be revised to take into account this new information. 
In this case the posterior probability found above becomes the 
prior probability in the right side of Eq. [10] and new posterior 
probabilities can be calculated. 

Here, we assume K - 5, representing the five groups of 
Peimbert I, Ila, lib. III and IV. At first only PNe that are pre- 
cisely classified in each of those groups are selected, in order to 
estimate the likelihood p(x\k, I) of the corresponding group. As 
a consequence, only nebulae that have all the four classificatory 
parameters are used. In the following discussion, we will leave 
the I dependence implicit in order to simplify the notation. 

Only seven objects from our sample were unambiguously 
classified as type IV planetary nebula. In order to improve the 
statistical analysis of type IV PNe, we have added to this group 
four other objects that were classified as IVC/IIbB. A similar 



problem occurs for type III PNe, which has only three unam- 
biguously classified objects. We added to this group 15 objects 
that were classified as IIIC/IIbB or IlbC/IIIB. This was done for 
the sole purpose of increasing the sample in order to not com- 
pletely hinder the statistical analysis for these less populated 
groups. However, by mixing PNe having uncertain classification 
with well-classified ones, the estimated multivariate likelihood 
for types III and IV PNe are surely less well-defined than those 
for the other types. 

In Table |3] we list the PNe used in the calculation of p{x\k) 
to each of the Peimbert types. In Table |4] we list the statisti- 
cal parameters of p(x\k). The mean of the variable x, is repre- 
sented by n(xi), and the notation c(xi, Xj) indicates the covari- 
ance between the variables x,- and xj. Accordingly, c{xi,Xi) is 
the variance of x/. Note that c(x,-,xy) = c{xj,Xi). Figure [1] shows 
two-dimensional projections of the four-dimensional likelihoods 
p{x\k) overimposed on the data for the PNe listed in Table [3] 
The ellipses mark the Icr confidence level for the bivariate dis- 
tribution of the parameters plotted, for each of the K groups. In 
this figure, the role of each classificatory parameter can be con- 
sidered. He/H and logN/0 are discriminant parameters for the 
earlier Peimbert types (I to II), but cannot separate between lib 
to IV PNe types. To separate these types, we need to consider 
|z| and |Ay|, which, on the other hand, are not good discrimi- 
nant parameters for the earlier types. We have verified that the 
Peimbert groups do not naturally arise as optimally decomposed 
groups of the four-dimensional variable mixing when a proper 
statistical method for the decomposition of population mixing is 
applied to the datc0. That is, the Peimbert classification is rather 
ad hoc, in the sense that it is based on the apriori expectation 
that objects within given parameter ranges are members of given 
stellar populations of the Galaxy. 

The shape of the ellipses in Figure[T]reveals some interesting 
properties of the Peimbert types. For instance, there is a trend of 
finding larger logN/O and |Ay| with increasing \z\ among type 
IV PNe. Because these PNe are thought to be members of the 
Galactic halo, we can conclude that this trend should be found 
among halo stars as well. The increase of |Ay| with increasing \z\ 
is indeed verified among halo stars. It indicates that stars with in- 
creasing height above the Galactic plane have a tendency to pro- 
gressively lag behind the Local Standard of Rest. On the other 
hand the logN/0-|z| trend found for these nebulae reinforces 
that metal-poor stellar populations could have a high N/0 abun- 
dance ratio, as discussed in section f32\ above. Other trends re- 
vealed by the ellipses are somewhat unexpected, as the decrease 
in |Ay| for increasing |z| among type III PNe. It is presently un- 
clear whether this trend is real or artificially imposed into the 
classification of type III PNe by the fuzzy limit in |z| for the sep- 
aration among type III to type IV PNe in the original Peimbert 
Classification. This trend deserves a further consideration after 
more data becomes available. 

The marginal likelihoods tell us about how each variable af- 
fects the classification of the nebulae. In Fig. |2]we show how 
the posterior probabilities pik\xq) depends on each variable. In 
those plots, the horizontal line correspond to p{k\xq) = 0.75, that 
we arbitrarily assumed as the limit for a safe classification of a 
nebulae in a given group. In a classification system that has few 
ambiguities, we expect to find almost each curve rising above 
p{k\xq) - 0.75 in a given variable range. Contrary to this expec- 
tation, Fig.|2]shows that the several p{k\xq) curves overlap each 
other considerably for some variables. As a consequence, only a 



- We have used the EMMIX software provided by McLachlan et al. 
(1999) 
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Type I 



BV5-1 


He2-15 


He2-143 


He2-153 


K3-70 


Ml-8 


Ml-51 


Ml-75 


NGC 2899 


NGC 6302 


NGC 6445 


NGC 6537 


NGC 6781 


NGC 6803 


NGC 7293 








Type Ila 


He2-29 


He2-37 


He2-51 


He2-86 


He2-119 


He2-140 


He2-141 


IC4I8 


IC 2149 


IC 2448 


IC 2501 


IC 2621 


K3-68 


Ml-4 


Ml-7 


Ml-13 


Ml-57 


Ml-58 


Ml -60 


Ml-61 


Ml -74 


Ml -79 


Ml-80 


M2-2 


M2-4 


M2-46 


M2-55 


M3-1 


M3-4 


M3-5 


M3-29 


MyCnl8 


Mzl 


Mz2 


NGC 650 


NGC 2022 


NGC 2346 


NGC 2371-72 


NGC 2438 


NGC 2440 


NGC 2452 


NGC 2792 


NGC 2867 


NGC 3132 


NGC 3195 


NGC 3699 


NGC 3918 


NGC 5189 


NGC 5315 


NGC 5882 


NGC 6572 


NGC 6720 


NGC 6741 


NGC 6818 


NGC 6881 


NGC 6884 


NGC 7026 


NGC 7027 


NGC 7354 


Th2-A 








Type lib 


H3-75 


Hbl2 


He2-115 


He2-47 


He2-55 


Hul-l 


Hu2-1 


IC 2003 


IC 2165 


IC 351 


IC4191 


IC 4406 


IC 4634 


IC 4776 


IC 5117 


J900 


K2-I 


M1-I4 


Ml-17 


Ml-5 


Ml-50 


Ml-6 


M3-54 


M3-6 


NGC 40 


NGC 3211 


NGC 3242 


NGC 4361 


NGC 6210 


NGC 6309 


NGC 6578 


NGC 6790 


NGC 6826 


NGC 6879 


NGC 6891 


NGC 6905 


NGC 7662 


Pel-1 


Tel 














Type III 


Fgl 


He2-5 


He2-21 


He2-67 


He2-99 


He2-149 


He2-157 


He2-158 


IC 4846 


J320 


KFL 19 


Ml-9 


M2-9 


M2-53 


M3-33 


NGC 1535 


NGC 6326 


PC 14 


Type IV 


BoBnl 


H4-1 


He2-I18 


IC 3568 


IC 4593 


K648 


PC12 


DdDml 


Me2-1 



Snl Vyl-2 



Table 4, Statistic parameters of p(x\k) 





Type I 


Type Ila 


Type lib 


Type III 


Type IV 


//(He/H) 


0.1503 


0.1135 


0.1035 


0.1004 


0.0983 


y"(logN/0) 


-0.0160 


-0.3973 


-0.8110 


-0.9161 


-0.7873 


y"(kl) (kpc) 


0.0812 


0.2251 


0.3761 


0.6859 


5.4504 


MIAl/|)(kms-') 


25.1667 


28.1200 


23.3615 


69.6389 


115.5273 


c(He/H,He/H) 


0.0003 


0.0003 


0.0002 


0.0003 


0.0001 


c(He/H, logN/O) 


0.0015 


-0.0005 


0.0001 


-0.0001 


0.0023 


c(He/H, Izl) 


-0.0002 


-0.0005 


-0.0000 


-0.0026 


0.0277 


c(He/H, |AV|) 


0.0168 


-0.0395 


-0.0146 


0.0642 


0.0143 


c(log N/0, logN/O) 


0.0444 


0.0373 


0.0300 


0.0275 


0.2231 


c(logN/0,k|) 


-0.0062 


-0.0071 


-0.0013 


-0.0331 


1.7941 


c(logN/0, |AV|) 


-0.7432 


0.0729 


-0.2483 


1.2787 


14.5046 


c(\z\,\z\) 


0.0034 


0.0445 


0.0663 


0.2045 


25.8489 


c(|z|,|AV|) 


-0.0391 


0.1964 


0.2333 


-5.9943 


218.0222 


c(\AV\,\AV\) 


203.3552 


282.4165 


269.3030 


495.8237 


8091.3482 


p(k) 


0.105 


0.419 


0.273 


0.126 


0.077 



small number of curves do reach p(k\xg) > 0.75. This indicates 
that Peimbert classification is rather ambiguous if only one clas- 
sificatory parameter is known. However, the same problem af- 
fects other population classification schemes in Astronomy, like 
the Galactic stellar population classification (see, for instance, 
Nemec & Nemec 1993). 

Seven intervals can be defined from these curves, within 
which the PNe can be classified with reasonable accuracy from a 
single parameter. Those are summarized in Table|5] We conclude 
that log N/0 is an important determinant parameter in the classi- 
fication of the types I and Ila, with two intervals of unambiguous 
classification in Table |5] |z| and |Ay| are also discriminatory pa- 
rameters between Type III and IV PNe, but not for those younger 
nebulae associated with the thin disk. He/H is not an unimportant 



parameter, but the overlap between the several p(k\}ie/H) curves 
hinders a secure nebulae classification unless He/H > 0.147. 

More narrow intervals can be defined by considering that 
only two or three classificatory parameters are known. We do 
not show the corresponding plots in this case, since no new in- 
formation can be taken from them that is not already present in 
the Figs. [T] and |2] If more than one classificatory parameter is 
known, it is best to use Eq.[TO]to fully classify the nebula. 

The accuracy of the likelihoods p(x\k) used in the posterior 
classifications depends on the accuracy of the prior classification 
that defined the 'fiducial' PNe samples. We had to incorporate 
some ambiguously classified Type III and IV PNe to these sam- 
ples in order to have a meaningfully large sample. On account 
of this, some ambiguity is likely to have been incorporated into 



10 C. Quireza et al.: A reanalysis of the Peimbert planetary nebulae types 



■ type I O type I la X type Mb s type III ► type IV 




0.06 0.08 O10 012 014 016 018 -1.50 -1.25 -1.00 -0.75 -O50 -025 0.00 025 20 40 60 80 100 120 140 160 180 SDO 220 240 



He/H log (N/O) |AV| (km/s) 

Fig. 1. Projections of the multivariate distribution of classificatory parameters in the Peimbert classification scheme. Different sym- 
bols are used for PNe unambiguously classified in the five main Peimbert group, as shown in the caption. An ellipse marks the 
1-sigma confidence level of each group bivariate distribution in the four-dimensional space (He/H, log N/O, |Ay|). While some 
groups are reasonably well separated in some of the plots, particularly He/H x log N/O and log N/O x |Ay|, they overlap consider- 
ably in other plots, on account of the ambiguity in the Peimbert classification and of real overlaps of some of the parameters that 
characterize the different nebular populations. 



Table 5. Best parameter intervals for safe classification of PNe 



Interval 


Class 


He/H > 0.147 


type 1 PNe 


log N/O > -0.005 


type 1 PNe 


-0.513 < log N/O < -0.324 


type 11a PNe 


0.95 < Izl < 1.60 kpc 


type 111 PNe 


Izl > 1.99 kpc 


type IV PNe 


66 < |AV| < 101 kms^' 


type 111 PNe 


|AV| > 121 kms^' 


type IV PNe 



p(x\k). A glimpse of Fig.[T]shows that the present 'fiducial' Type 
III PNe sample comprises PNe with lower He/H abundance than 
the 'fiducial' Type IV PNe sample, which are supposed to be 
older than Type III PNe and, thus, more He poor This problem 
can be eliminated once the number of prior classified nebulae 
increases as more data becomes available^ 

3 Recently, Pereira & Miranda (2007) identified PNG 232.0-t-05.7 as 
a halo object, increasing the small group of known halo planetary neb- 
ulae (PNG 232.0-f05.7 present a very low N/O abundance ratio, with 



In Table|6](available also in electronic form) we list the pos- 
terior probability p(k\y) for the membership of each planetary 
to the Peimbert groups, where y is the vector composed by the 
classificatory parameters of the nebula that are known. For com- 
parison purpose, we also list the pre-classification, as in Table[T] 
While the posterior classification agrees reasonably well with 
the prior classification for those PNe used in the fiducial' sam- 
ples, the novelty of this result is that we can now give a quanti- 
tative estimate of how likely each ambiguously classified nebula 
belongs to a given Peimbert group. This result is particularly im- 
portant in the case of nebulae for which a few classificatory pa- 
rameters are known. The resulting posterior classification listed 
in the eighth column corresponds to the group with higher pos- 
terior probability. Compared to the percentages provided in the 
§ l3.1l for the numbers of PNe in each group, we have now that, 
excluding the objects that were pre-classified as bulge PNe (28% 
of the whole sample), 42% of the whole sample of 476 PNe are 
type II PNe (31% Ila and 11% lib), 11% are type III, 9% are 



log(N/0) = -0.89). Analogously, PNG 034.5-11.7 was identified as a 
type 111 planetary nebula. Unfortunately none of these new identified 
PNe have He/H measurements. 
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|z| (kpc) |aV| (km/s) 

Fig. 2. Posterior probability, p(k\xg), or the fraction of PNe which are likely to be members part of a group k, when only one single 
information for the nebula is known. The horizontal line indicates pik\xg) = 0.75. According to the figure, log N/0 and |z| are the 
most important parameters in the Peimbert classification. 

type I, and 6% are type IV. There are still a number of PNe with 
less than 50% of chance to pertain to a certain group, including 3 
objects that still cannot be classified, having equal probabilities 
to belong to two different types. Those PNe represent 4% of the 
sample. 

In order to avoid indiscriminate use of the results presented 
in Table |6] we wish to present some critical comments. First of 
all, we find that a number of 30 PNe are now classified as halo 
type IV. Those include objects classified as halo PNe in the lit- 
erature (Howard et al. 1997) plus objects with ambiguous pre- 
classification among two or more Peimbert types. We have no- 
ticed, however, that this group contains 1 1 PNe previously clas- 
sified as type Ila. Although the scale heights for most of these 
objects are less than 1 kpc, their peculiar velocities are high. The 
inclusion of these objects in the type IV group may be partially 
explained by the large scatter in N/O observed in halo objects. 
Most of the pre-classified Ila PNe have low He/H (< 0.125) but 
log(N/0) > -0.60. 

It also should be noted that in Table |6] we list posterior clas- 
sification for type V bulge PNe. This was done in order to inves- 
tigate the main characteristics of this group, given the meaning 
of the sequence of types I to IV as being a sequence of intervals 
of progenitor mass (see discussion ahead). We do not recom- 
mend the use of our resulting posterior classification for those 
PNe pre-classified as type V for application in studies dealing 
with Galactic evolution or related purposes. In those case, bulge 
objects should be used as a separate group. 



Planetary nebula bulge population: It is interesting to inves- 
tigate the characteristics of the PNe in the bulge population us- 
ing the posterior classification of PNe that have a type V prior 
classification. Until now, we have discarded from the statistical 
analysis objects belonging to this class. This procedure has as- 
sured a better characterization of the types I to IV. Here we can 
study this class separately. We calculated the posterior classifi- 
cation of all PNe with prior classification VA, VD and VF in 
Table [1] Those are listed in Table |6] For this group of objects 
we have 6 type I, 41 type Ila, 19 type lib, 35 type III and 32 
type IV PNe. We conclude that type V PNe are, indeed, a mixed 
group comprising a variety of chemical composition and stellar 
masses. They are likely to be misclassified in all other Peimbert 
groups if the proper classificatory parameters for type V PNe 
are not taken into account. The frequency of misclassified type 
V PNe as Type III or IV PNe, however, compared to the prior 
probabilities p{k) is larger by a factor of at least 2. This suggests 
that most type V PNe are more similar to type III and IV PNe, in 
agreement with what we would expect from Galactic population 
studies, according to which the bulge population is closer to the 
thick disk and halo population than to the thin disk population 
(Feltzing & Gilmore 2000; FeiTeras et al. 2003). 

5. Discussion and Application 

In the Figs [T] and |2] we note significant overlap among the dis- 
tributions of the classificatory parameters, which best define the 
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different types I to IV. As commented in the previous sections, 
the lack of sharp boundaries between these intervals is one of 
the main causes for uncertainties in the classification of PNe 
into these types. The overlap in the pix\k) curves should be re- 
duced with a better characterization of the sample, which can 
be reached with the improvement of the observational data, tak- 
ing into account new and more precise measurements. A larger 
number of PNe pre-classified into the types III and IV is also 
of major importance, since those groups are poorly populated 
in our analysis. We do not expect, however, that the overlaps in 
Figs. [1] and |2] will be completely removed. Some overlaps are 
probably real, since different nebula populations have common 
properties. Some discussion was already presented in § |3.2| con- 
cerning the large scatter observed in abundances from thin disk 
to halo, and possible effects of the motion of the progenitor stars 
in the position of the PNe. Thus, overlaps should exist not only in 
the curves of the He/H and log N/0 abundance ratios, but also, 
although in a smaller range of coincidence, in the curves for |Ay | 
and \z\, particularly for the older PNe. 

This discussion reinforces our idea that statistical methods 
can greatly contribute to a better characterization of groups of 
PNe. It is also clear that these groups cannot be well defined 
using a limited number of criteria, and the consideration of fur- 
ther quantities that can be related to stellar evolution, kinematics 
or spatial distribution can probably help in the characterization 
of stellar populations. Following this line of thought, Phillips 
(2005) has evaluated trends in several planetary nebula proper- 
ties for the Peimbert types I, Ila, lib, and III, finding marked 
differences in properties among the four Peimbert classes, and 
confirming a likely differentiation in progenitor and ionized shell 
masses. Among the quantities analyzed were the: nebular radii, 
expansion velocity, density, shell ionized masses, morphology, 
radial velocity. Galactic latitude, 5 GHz brightness temperature, 
Zanstra temperature (likely to be related to central star effec- 
tive temperature), central star luminosity, He^^/He^, dust-to-gas 
mass ratios, and molecular-to-ionized mass ratios. In the follow- 
ing subsections we present additional comments on the relation 
between morphological classes and PNe Peimbert types. We also 
briefly comment on isotopic ratios as indicators of the mass of 
planetary nebula progenitor star. 



5. 1 . Planetary nebula morphological class as a stellar 
population indicator 

Several studies use planetary nebula morphology as indicator of 
stellar populations. It was first noticed by Peimbert & Torres- 
Peimbert (1983) that a large fraction of type I PNe (He and N 
rich) show bipolar structure (see also Calvet & Peimbert 1983). 
Later on, the association of progenitors of bipolar PNe and the 
higher-end mass range of the AGB stars was confirmed in a se- 
ries of papers dealing with the correlation of the planetary neb- 
ula morphology and properties of the nebula and its central star 
(Amnuel 1995; Corradi et al. 1997; Corradi & Schwarz 1995; 
Gornyetal. 1997; Stanghellini et al. 1993, 2002). Morphological 
classes are usually defined on the basis of the Ha, [O m] A 5007 
or [Nn] /1 6584 images. Attempts to establish a morphological 
classification have been done by many authors (Greig 1972; 
Zuckerman & Afler 1986; Balick 1987; Chu et al. 1987; Schwai'z 
et al. 1992; Manchado et al. 1996). Some variation in the def- 
inition of groups and subgroups exist, depending on the level 
of detail which is taken into account, by considering or not ad- 
ditional structures and morphological features (inner structures, 
ansae, rings, multiple shell structures, and other). Also, different 
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Fig. 3. Distribution of morphological classes and PNe types (ac- 
cording to our posterior classification) for: a) 77 objects in com- 
mon with Stanghellini et al. (2002) and b) 77 objects in common 
with Gorny et al. (1997) and Stasiilska et al. (1997). The scales 
in the captions are normalized numbers: in the upper plot (a), 
groups were normalized by 10 type I, 39 Ila, 14 lib, 6 III and 8 
IV PNe; while in the lower plot (b) we have 9 type I, 3 1 Ila, 20 
lib, 13 III and 4 IV PNe. 



interpretation of a given image, may result in differences in the 
classification assigned by different authors. 

Fig. [3^ shows a comparison between the Peimbert types and 
the morphological classes in Stanghellini et al. (2002) for PNe in 
common with their study. Fig.[3l5 is analogous to Fig. [3^, but now 
showing PNe in common with Gorny et al. (1997) and Stasiilska 
et al. (1997). Peimbert types are given according to our posterior 
classification. The number of objects in each interval are plotted. 
For a proper comparative analysis, the groups are normalized 
by the total number of PNe in each of the Peimbert types. The 
gray color scale refers to the normalized numbers, where darker 
colors indicate higher numbers of objects. In Fig. [3^ we have 10 
type I PNe, 39 Ila, 14 lib, 6 III and 8 IV. In Fig. ^ we have 9 
type I, 31 Ila, 20 lib, 13 III and 4 IV, summing a total number of 
77 objects in each figure. 

In Gorny et al. (1997) and Stasiilska et al. (1997), PNe were 
classified as point-symmetric (P), bipolar (B), elliptical (E) and 
irregular (Ir). Further subdivisions were given to these groups, 
with the subdivision of the point-symmetric into PE and P, ac- 
cording to whether the knots are embedded or not in an overall 
structure. Analogous classification was applied to the bipolar, 
resulting in the BE and B sub-groups. The elliptical were subdi- 
vided into e, elliptical disk without structure; E, clear elliptical 
morphology; and EH, almost circular with very regular shape 
and bright rim (we have none of this last morphological class). 
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Stanghellini et al. (2002) used a more simplified morpho- 
logical scheme, taking into account three main morphological 
classes, namely round (R), elliptical (E) and bipolar (B), each 
morphological classes being associated to different stellar pop- 
ulations. Some details were included by considering the bipo- 
lar core group (BC). According to this classification, round PNe 
should be the remnants of the lower mass progenitors, ellipti- 
cal planetary nebulae represent low to intermediate mass pro- 
genitors, and bipolar PNe are the result of high stellar mass 
progenitor evolution (Stanghellini et al. 2002, 2006; Manchado 
2004). Hereafter, we will refer to this classification as Simplified 
Morphological Classification. 

From Fig. [3^, most of the type I coiTespond to bipolar ob- 
jects, most of the types Ila, lib and III, are elliptical and most of 
the type IV are round PNe, which roughly confirms the expected 
correlation between the morphological classes and the stellar 
population. A similar conclusion can be drawn from the appli- 
cation of a contingency table test for the null hypothesis that the 
Peimbert and the Simplified Morphological Classification are in- 
dependent, using the data from Fig. |3^. We have verified that 
this null hypothesis can be rejected with a 0.05 level of signif- 
icance (P-value = 0.03), i.e., the PNe classes in the Simplified 
Morphological Classification are not entirely uncorrected with 
those classes from the Peimbert Classification. However, there is 
no straight coiTespondence between these morphological classes 
and the Peimbert types. For example, for 17 nebulae in common 
with the round morphological group, 6 are Ila, which is not con- 
sistent with the supposition that round PNe were predominantly 
in the low-mass domain. 

In Fig. [3j), for the Gorny et al. (1997) plus Stasinska et al. 
(1997) sample, the interpretation of correspondence between the 
classifications is not clear, probably due to the many subdivi- 
sions in the morphological classes. Some speculation, however, 
can be done about this plot. Point-symmetric PNe are supposed 
to be present in a variety of main-body morphologies, which 
implies that in the Simplified Morphological Classification, 
they could be found in the round, bipolar or elliptical classes 
(StanghelHni et al. 2002). StanghelHni et al. (2002) call atten- 
tion to a similarity in the definition of bipolar core (BC) and 
the one used by Gorny et al. (1997) for embedded bipolar (BE), 
which, also, would be classified as R or E in the Simplified 
Morphological Classification. Roughly, we find that Type I ob- 
jects are almost equally distributed in the E and B classification. 
For all the other Peimbert types, the E morphology dominates. 
Since Gorny et al. (1997) do not distinguish between round and 
elliptical, it is possible that some of the type IV PNe could be 
classified as round. Thus, although less obviously, the distribu- 
tions in Fig. [3J3 still can be interpreted as showing a relation of 
morphological class and PNe mass progenitor, but the variety of 
morphological structures found inside of each of the types I to 
IV avoid strong statements. A contingency table test for the inde- 
pendence between the Gorny et al. morphological classification 
and the Peimbert Classification reinforces our conclusions: The 
f -value for rejecting the null hypothesis in this test is 0.52, indi- 
cating that both classification schemes have high probability of 
being independent. 

A N/O X He/H plot of elliptical, round and bipolar PNe 
would show significant overlap between different morphological 
samples. One example can be found in Stanghellini et al. (2006, 
Fig. 1), which shows elliptical and round PNe in the same data 
range, including some in the typical data range for bipolar PNe. 
Thus, although we do find indication that the main morphologi- 
cal classes, namely bipolar, elliptical and round are in some way 
coiTelated to the Peimbert types, which are supposed to corre- 
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Fig. 4. The '^C/'^C isotopic ratio in the different Peimbert types, 
according our post-classification. Regardless uncertainties in ob- 
servational data and theoretical models, large scatters are possi- 
ble for the isotopic ratio in each type. 



spond to a sequence of stellar populations of the Galaxy, it is 
clear that PNe morphological class cannot be used as unique in- 
dicator of stellar population. 

5.2. The '^C/'^C isotopic ratio as a stellar mass indicator 

The abundance content of several chemical elements and iso- 
topic ratios (^He, "^He, ^^C/^^C, '^N/'^N, '^O/'^O) in PNe re- 
sults from complex combination of different dredge-ups, mix- 
ing, mass loss and nuclear processing episodes that occur along 
the evolution of their progenitors. In particular, abundances that 
will be significantly affected should be those of He, C, and N 
(and possibly O). Those changes will depend mainly on the ini- 
tial stellar mass and, in a minor degree, on the metallicity. This 
leads to the possibility of using observations of isotopic abun- 
dance ratios in PNe as constraints of stellar mass and age. 

Standard models predict that in specific phases of the stel- 
lar evolution, the external convective envelope moves down- 
ward reaching internal layers where H burning has occurred, and 
bring up processed material to the stellar surface. As a result, 
changes in the surface abundance will occur after each dredge- 
up episode. A significant increase in the surface '^C/'^C isotopic 
ratio is expected after the third dredge up in an AGB, where the 
H burning shell is expected to be enriched with '^C due to a pre- 
vious He shell flash. The third dredge up results in high '^C/'^C 
ratios, larger than 20 (Boothroyd & Sackmann 1999; Renzini 
& Voli 1981). As observed values of the '^C/'-'C isotopic ra- 
tio on the RGB and later evolution (Charbonnel & Nascimento 
1998) suggest that values predicted by standard models are over- 
estimated, additional processes, employing deep extra mixing 
mechanisms during RGB and/or AGB phase are usually applied 
in order to explain discrepancies (Charbonnel 1994, 1995; van 
den Hoek & Groenewegen 1997; Marigo 2001). Those non- 
standard processes result in a significant enhancement of the '^C 
abundance in the surface layers. 

In low mass stars (^ 2 Mq), nonstandard mixing models 
that incorporate cool bottom processing (during RGB and/or 
AGB phase) are usually applied to reproduce observational data. 
In this process the bottom of the convective envelope remains 
cool while an ad hoc mixing mechanism brings material down 
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to layers hot enough for some nuclear processing. Such mod- 
els can lead the ^^C/'-'C isotopic ratio to values as low as 4-5 
(Wasserburg et al. 1995). For masses larger than x 4 Mq, stars 
develop hot bottom burning in the AGB phase, where deep con- 
vective envelop with very high base temperature activates the 
CN cycle, as it penetrates into the H burning shell, processing 
'^C into '^C and ''*N. This process tend to drive '^C/'^C toward 
~ 3.3 (Smith & Lambert 1990; Frost et al. 1998). 

In Fig.m we present the posterior classification of PNe with 
observed '^C/'-'C isotopic ratios. Data were taken from Balser et 
al. (2002), supplemented by Palla et al. (2000, 2002) and Rubin 
et al. (2004). For Palla et al. (2000) and Balser et al. (2002) mea- 
surements were made from millimeter transitions of '^CO and 
'^CO. These samples comprise a larger number of type I and 
mostly type 11a PNe. Later types, comprehending PNe originated 
from older stars, are found in larger (although still small) num- 
ber in the studies from Palla et al. (2002) and Rubin et al. (2004), 
which used the C iii] multiplet, containing the '^C lines at 1906.7 
and 1908.7 A and the very weak '^^C fine at 1909.6 A. For a 
large number of PNe, in which the '^'C line was not detected, 
only lower limits of the '^C/'^C isotopic abundance ratios were 
provided. Most of these PNe are elliptical or bipolar (one bipolar 
core), regardless of the origin of the sample. 

The types I to IV should be interpreted as a sequence of de- 
creasing mass for the PNe progenitor star in the main sequence. 
Those estimates should provide a critical check of the Peimbert 
classification as a sequence of mass intervals. It is beyond the 
scope of this paper to analyze these types as a function of stellar 
mass, since there are large uncertainties associated with those 
estimates, and no interesting constraint (no new information) 
would be provided, in addition to the ones already discussed in 
Balser et al. (2002) and referred works that used mass estimates. 
Some attempts have already been made, in those studies, to es- 
tablish a relationship of the '^C/'^^C isotopic ratio in PNe with 
the mass of the progenitor star in the main sequence, but large 
uncertainties in the determination of the stellar mass and incon- 
sistencies in standard stellar evolution theory and observations 
hinder any strong statements (Balser et al. 2002). The sequence 
of types in Fig.|4]should provide an independent qualitative pic- 
ture of any possible trend of the '^C/'^C ratio with the mass of 
the progenitor star. 

In the figure, the observations of the CO millimeter-wave and 
Cm] transition in PNe yield '^C/'^C values between ~ 2.2 and 
31, with progenitor masses from ~ 1 to 4 Mq. Lower limit data 
expand these intervals to ~ 0.5 and 45, roughly, for the same 
mass interval. Low values (lower than 5) are observed for the 
'^C/'^C isotopic ratios in the older types III and IV, in agree- 
ment with results provided by models, that use deep mixing cool 
bottom processing. Because a large number of lib PNe may have 
masses lower than 2 Mq, cool bottom burning should also occurs 
in, at least, part of the group of PNe classified as type lib. In the 
intermediate mass interval, comprising types Ila and lib, a large 
scatter in the '^C/'^'C ratio is observed, suggesting that only part 
of the stars suff'er deep mixing. In fact, a problem remains for 
stars between 2-4 Mq, since cool bottom processing, is only ex- 
pected to occur in progenitor stars of < 2 Mq. Thus, low '^C/'"^C 
between 2 and 4 Mq might indicate further processing on the 
AGB or that the observed ratios are underestimated. The type 
I planetary nebula in the figure shows a low '^C/^^C ratio, in 
agreement with models that use hot bottom burning. This ex- 
tra mixing mechanism, however, should occur only in the more 
massive (<: 4 Mq) objects in the type I mass range. Because 
the type I PNe group also contain stars with masses lower than 



4 Mq, we should also observe higher values for the '^C/'^C ratio 
within this group. 

In summary, we could expect to observe a trend of increas- 
ing '^C/'^C ratio in the direction of the lower masses type I PNe, 
with some scatter A large scatter should be observed within the 
types Ila and lib, with a slight trend of decreasing '^C/'^C ra- 
tios in the direction of the lower mass lib PNe. Because some 
scatter in the isotopic ratio should always be observed in the 
types I to lib, the only significant constraint the '^C/'^C isotopic 
ratio could provide would be for the types III and IV, which 
are expected to present low values. These conclusions, how- 
ever, cannot be taken for granted, since larger uncertainties are 
still associated with stellar evolution models and observational 
data. Observations of NGC 3242 (estimated mass of 1 .2 Mq) 
by Palla et al. (2002) resulted in a carbon isotopic ratio > 38, 
in agreement with standard stellar models. This could provide 
evidences that not all progenitors undergo a phase of deep mix- 
ing during the last stages of its evolution. Based on chemical 
evolution models, more than 90% of low-mass stars must un- 
dergo cool bottom processing (Charbonnel et al. 1998). Thus, 
NGC 3242 must belong to the small group of standard low-mass 
stars, consisting of ~ 10% of stars with mass lower than 2 Mq. 
Additionally, a large fraction of the data presented is lower lim- 
its. A large number of observational measurements, in particu- 
lar from C III] transition, in order to cover older PNe, is desir- 
able. Concerning theory, new models that incorporate the effects 
of axial rotation (Meynet & Maeder 2002; Meynet et al. 2006) 
or of variable molecular opacities (opacities consistent with the 
changing photospheric chemical composition) in the competi- 
tion between third dredge-up and hot bottom burning during 
thermally pulsating AGB phase (Marigo 2007; Marigo & Girardi 
2007) can provide progressive advances in the present under- 
standing of stellar evolution in the intermediate mass range. 



6. Summary 

We have re-analyzed the criteria used to characterize the 
Peimbert types I, Ila, lib. III and IV, by performing a Bayesian 
statistical classification of a large sample of PNe previously clas- 
sified into these groups. Our sample consists of 476 PNe, whose 
nebular properties have been compiled from selected sources 
from the literature. A subsample, composed only by PNe satis- 
factorily distributed into the types I to IV, was used to calculate 
the joint probability density function of the classificatory param- 
eters (He/H, log N/O, |z| and |Ay|) in the five groups considered. 
These distributions were then used to calculate the probability 
of each planetary nebula to be member of a given Peimbert type. 
This probability, called posterior probability, provides a quanti- 
tative result of the representativity of the object within its group. 

Our posterior classification increases the number of PNe 
classified into the Peimbert types, extending the classification 
to PNe that are ambiguously classified in the traditional method. 
Even PNe for which a few classificatory parameters are known 
can be classified. 

Uncertainties in our posterior classification are given by the 
posterior probabilities for an object to be member of each of the 
types I, Ila, lib. III and IV. The classification can be improved 
as new observational data is added to the sample, and the clas- 
sification of the nebula is revised, to take into account this new 
information. 

Our analysis confirms that: i) The types I and Ila are well de- 
fined in terms of abundance, contrary to the other three groups 
lib. III and IV. In particular, log N/O is more relevant to define 
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the different types, ii) All the types are relatively well character- 
ized by the height, which clearly is one of the most impor- 
tant parameters to classify a halo PNe. iii) The peculiar velocity 
tends to increase as we go from thin disk, thick disk to halo, as 
expected, iv) \z\ and |Ay| are essential parameters in the classifi- 
cation of the different components of the Galaxy (thin disk, thick 
disk and halo), although they cannot be used alone to distinguish 
the three types pertaining to the thin disk (I, Ila and lib). 

We found that one of the main causes for ambiguities exist- 
ing in the original method of selection of PNe into the types I to 
IV should be associated with the difficulty in defining sharp mass 
boundaries for each of these groups. This conclusion can be in- 
ferred from the superposition we have found of the curves of 
distributions of the classificatory parameters (He/H, log N/0, \z\ 
and I Ay I), which best define the different types I to IV. The over- 
lap in the probability density functions will be reduced as the 
observational data is improved by taking into account new and 
more precise measurements (not available for all objects) and in- 
cluding a larger number of PNe pre-classified into the types III 
and IV, since those groups are poorly populated in our analysis. 

The overlaps in the curves of distributions of He/H, log N/O, 
\z\ and I Ay I, however, should be in part explained by real overlap 
of some properties of different stellar populations. This suggests 
that the classification can be improved if a larger number of clas- 
sificatory parameters is taken into account. 

We have checked how morphological PNe classifications 
compare to the Peimbert Classification. We have found that the 
SimpUfied Morphological Classification used in StanghelUni et 
al. (2006) is not uncorrected with the Peimbert Classification. 
There is a trend of having more bipolar type I PNe, elliptical 
type n and HI PNe and round type IV planetary nebula in the 
Galaxy. This indicates that PNe morphology may be an indica- 
tor of the stellar population to which the planetary nebula be- 
longs. However, the more detailed morphological classification 
by Gorny et al. (1997) has been found to be independent from the 
chemokinematical Peimbert Classification, probably on account 
of the larger variety of morphologies considered. 

Finally, we have used the posterior PNe classification devel- 
oped in this paper to explore whether there is a trend between the 
'^C/'-'C ratio and the mass of the progenitor star, which is likely 
to be correlated with the Peimbert types if these correspond to a 
sequence of stellar populations of the Galaxy. Our results do not 
rule out the existence of this trend, although the available data 
still prevent us from drawing strong conclusions. 
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Table 1. Parameters and pre-classiflcation (full table and reference lists are available electronically) 
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* Radial velocities from the catalog of Durand et al. (1998). For every object, the heliocentric radial velocities have been converted to the Local Standard of Rest radial velocities 
(see text). 

References of abundances. (3) Costa et al. 2004; (21) Perinotto 1991; (22) Perinotto et al. 1994; . . . 
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References of the 5 GHz flux densities. (9) Zhang 1995 (see references in); . . . 



►13 
I 



ft 



o 
s. 

re 

I 



a 

So 



C. Quiieza et al.: A reanalysis of the Peimbert planetary nebulae types, Online Material p 3 



Table 6. Posterior probabilities for the planetary nebulae to belong to Peimbert classes 
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0.00 


1.00 


0.00 


0.00 


0.00 


VF 


Ila 


Hl-29 


0.01 


0.59 


0.28 


0.02 


0.10 


(V/IIa)D 


Ila 


Hl-30 


0.00 


1.00 


0.00 


0.00 


0.00 


VF 


Ila 


Hl-31 


0.00 


0.21 


0.05 


0.73 


0.01 


VF 


III 


Hl-32 


0.00 


0.01 


0.11 


0.87 


0.01 


VF 


III 


Hl-33 


0.00 


0.98 


0.02 


0.00 


0.00 


VF 


Ila 


Hl-34 


0.00 


1.00 


0.00 


0.00 


0.00 


VF 


Ila 


Hl-35 


0.00 


0.00 


0.00 


0.23 


0.77 


VF 


IV 


Hl-39 


0.00 


0.00 


0.00 


0.00 


1.00 


VF 


IV 


HI -40 


0.00 


0.82 


0.02 


0.15 


0.01 


VF 


Ila 


Hl-41 


0.00 


0.00 


0.00 


1.00 


0.00 


VF 


III 


Hl-42 


0.00 


0.96 


0.04 


0.00 


0.00 


VF 


Ila 


Hl-43 


0.00 


0.21 


0.00 


0.79 


0.00 


VF 


III 


HI -44 


0.00 


0.22 


0.00 


0.34 


0.44 


VD 


IV 


HI -45 


0.00 


0.00 


0.00 


0.00 


1.00 


VF 


IV 


HI -46 


0.00 


0.00 


0.00 


0.46 


0.54 


VF 


IV 


Hl-47 


0.00 


0.13 


0.00 


0.87 


0.00 


VF 


in 


Hl-50 


0.00 


0.00 


0.00 


0.99 


0.01 


VF 


m 


Hl-51 


0.00 


0.00 


0.00 


0.00 


1.00 


VD 


IV 


Hl-54 


0.00 


0.00 


0.00 


0.03 


0.97 


VF 


IV 


Hl-55 


0.01 


0.42 


0.42 


0.10 


0.05 


indef. 


Ila/IIb 


Hl-56 


0.00 


0.00 


0.00 


0.83 


0.17 


VF 


III 


Hl-58 


0.00 


0.00 


0.10 


0.90 


0.00 


(V/nb/lll/lV)D 


III 


Hl-59 


0.00 


0.01 


0.99 


0.00 


0.00 


VF 


lib 


Hl-60 


0.00 


0.41 


0.59 


0.00 


0.00 


VF 


lib 


Hl-63 


0.00 


0.00 


0.02 


0.98 


0.00 


VF 


m 



C. Quiieza et al.: A reanalysis of the Peimbert planetary nebulae types, Online Material p 4 



Table 6. continued. 



Name 


/j(type I[y) 


/7(type na[y) 


/>(type nbly) 


/)(type m\y) 


p(type IVLy) 


Pre-class 


Post-class 


Hl-66 


0.00 


1.00 


0.00 


0.00 


0.00 


VF 


Ila 


HI -67 


0.00 


0.00 


0.00 


0.00 


1.00 


VF 


IV 


H2-1 


0.00 


0.01 


0.99 


0.00 


0.00 


VF 


lib 


H2-10 


0.00 


0.00 


0.00 


0.00 


1.00 


VF 


IV 


H2-11 


0.00 


1.00 


0.00 


0.00 


0.00 


VF 


Ila 


H2-15 


0.00 


0.00 


0.00 


0.00 


1.00 


VF 


IV 


H2-18 


0.00 


0.15 


0.00 


0.00 


0.85 


VF 


IV 


H2-24 


0.00 


0.00 


0.00 


0.00 


1.00 


VF 


IV 


H2-25 


0.00 


0.11 


0.89 


0.00 


0.00 


VF 


lib 


H2-26 


0.00 


0.00 


0.00 


1.00 


0.00 


VD 


III 


H2-37 


0.00 


0.00 


0.00 


0.75 


0.25 


VD 


III 


H2-46 


0.07 


0.87 


0.05 


0.00 


0.01 


VE/IIaD 


Ila 


H2-48 


0.00 


0.00 


0.00 


1.00 


0.00 


IIbC/(IlI/IV)B 


III 


H3-75 


0.00 


0.08 


0.92 


0.00 


0.00 


IlbF 


lib 


H4-1 


0.00 


0.00 


0.00 


0.00 


1.00 


IVF 


IV 


Hb4 


0.02 


0.96 


0.02 


0.00 


0.00 


VA 


Ha 


Hb5 


0.99 


0.01 


0.00 


0.00 


0.00 


VA 


I 


Hb6 


0.04 


0.95 


0.01 


0.00 


0.00 


VA 


Ila 


Hb8 


0.00 


0.00 


0.00 


1.00 


0.00 


IlaE 


III 


Hbl2 


0.00 


0.07 


0.93 


0.00 


0.00 


nbF 


nb 


He2-5 


0.00 


0.00 


0.00 


0.72 


0.28 


niF 


m 


He2-7 


0.00 


0.49 


0.01 


0.49 


0.01 


IlaC 


Ila/III 


He2-9 


0.27 


0.61 


0.10 


0.02 


0.00 


indef. 


Ila 


He2-15 


1.00 


0.00 


0.00 


0.00 


0.00 


IF 


I 


He2-21 


0.00 


0.00 


0.00 


1.00 


0.00 


IlbC/IIIB 


III 


He2-28 


0.00 


0.75 


0.13 


0.12 


0.00 


(na/IIb)D* 


Ila 


He2-29 


0.02 


0.92 


0.06 


0.00 


0.00 


IlaF 


Ila 


He2-37 


0.03 


0.80 


0.17 


0.00 


0.00 


IlaF 


Ila 


He2-47 


0.00 


0.21 


0.79 


0.00 


0.00 


IlbF 


lib 


He2-48 


0.00 


0.74 


0.24 


0.02 


0.00 


IlaD 


Ila 


He2-51 


0.00 


0.99 


0.01 


0.00 


0.00 


IlaF 


Ila 


He2-55 


0.00 


0.43 


0.57 


0.00 


0.00 


IlbF 


lib 


He2-57 


0.00 


0.56 


0.38 


0.06 


0.00 


IlaD 


Ila 


He2-67 


0.00 


0.00 


0.00 


1.00 


0.00 


IlbC/IIIB 


III 


He2-71 


0.00 


0.00 


0.00 


1.00 


0.00 


IlbC/IIIB 


III 


He2-76 


0.00 


0.35 


0.38 


0.27 


0.00 


IlbD 


lib 


He2-77 


0.00 


0.00 


0.90 


0.10 


0.00 


IlbD 


lib 


He2-86 


0.02 


0.98 


0.00 


0.00 


0.00 


IlaF 


Ila 


He2-97 


0.00 


0.00 


0.00 


0.46 


0.54 


nbc/(in/iv)B 


IV 


He2-99 


0.00 


0.00 


0.02 


0.97 


0.01 


IlbC/IIIB 


III 


He2-105 


0.00 


0.00 


0.46 


0.53 


0.01 


IlbD 


III 


He2-108 


0.00 


0.00 


0.00 


0.96 


0.04 


indef. 


III 


He2-lll 


1.00 


0.00 


0.00 


0.00 


0.00 


IC 


I 


He2-112 


0.00 


1.00 


0.00 


0.00 


0.00 


ID 


Ila 


He2-115 


0.00 


0.11 


0.89 


0.00 


0.00 


IlbF 


lib 


He2-117 


0.42 


0.57 


0.00 


0.01 


0.00 


IC 


Ila 


He2-118 


0.00 


0.00 


0.00 


0.77 


0.23 


IVC/IIbB 


III 


He2-119 


0.13 


0.87 


0.00 


0.00 


0.00 


IlaF 


Ila 


He2-123 


0.36 


0.64 


0.00 


0.00 


0.00 


IC 


Ila 


He2-131 


0.00 


0.92 


0.03 


0.05 


0.00 


indef. 


Ila 


He2-137 


0.00 


0.95 


0.02 


0.03 


0.00 


IlaF 


Ila 


He2-138 


0.00 


0.87 


0.13 


0.00 


0.00 


indef. 


Ila 


He2-140 


0.00 


0.90 


0.10 


0.00 


0.00 


IlaF 


Ila 


TJaO 1/11 


u.u/ 


f\ CO 


U.U4 


f\ f\f\ 


f\ f\f\ 

u.uu 


liar 


Ila 


He2-143 


0.98 


0.02 


0.00 


0.00 


0.00 


IF 




He2-149 


0.00 


0.00 


0.00 


1.00 


0.00 


IlbC/IIIB 


III 


He2-152 


1.00 


0.00 


0.00 


0.00 


0.00 


IF 




He2-153 


0.97 


0.03 


0.00 


0.00 


0.00 


IF 




He2-157 


0.00 


0.01 


0.02 


0.96 


0.01 


IlbC/IIIB 


III 


He2-158 


0.00 


0.00 


0.45 


0.54 


0.01 


IIIC/IIbB 


III 
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Table 6. continued. 



Name 


/j(type I[y) 


p{type na[y) 


p(type nbly) 


/)(type m\y) 


p(type IVLy) 


Pre-class 


Post-class 


He2-164 


0.38 


0.44 


0.18 


0.00 


0.00 


indef. 


Ha 


He2-170 


0.00 


0.00 


0.00 


0.00 


1.00 


IIbC/(III/IV)B 


IV 


He2-175 


0.00 


0.00 


0.00 


0.18 


0.82 


IlaC 


IV 


He2-250 


0.00 


0.00 


0.00 


0.96 


0.04 


VF 


III 


He2-260 


0.00 


0.00 


0.00 


1.00 


0.00 


IIIC/IIbB 


III 


He2-262 


0.00 


0.00 


0.00 


0.00 


1.00 


VF 


IV 


He2-406 


0.13 


0.84 


0.03 


0.00 


0.00 


(IlaA^)D 


Ila 


He2-434 


0.00 


0.09 


0.66 


0.24 


0.01 


(nb/III/IV)D 


lib 


He3-1357 


0.00 


0.06 


0.59 


0.22 


0.13 


(IIb/III/IV)D 


lib 


Hf2-2 


0.00 


1.00 


0.00 


0.00 


0.00 


VD 


Ila 


Hul-l 


0.00 


0.03 


0.96 


0.00 


0.01 


IlbF 


lib 


Hul-2 


0.00 


1.00 


0.00 


0.00 


0.00 


IB 


Ila 


Hu2-1 


0.00 


0.06 


0.93 


0.01 


0.00 


IlbF 


lib 


Hu4 


0.36 


0.64 


0.00 


0.00 


0.00 


(V/I)D 


Ila 


Hu5 


0.99 


0.01 


0.00 


0.00 


0.00 


(V/I)D 


I 


Hu6 


0.11 


0.81 


0.00 


0.00 


0.08 


(V/IIa)D 


Ila 


IC 351 


0.00 


0.04 


0.96 


0.00 


0.00 


IlbF 


lib 


IC418 


0.00 


0.62 


0.32 


0.06 


0.00 


IlaF 


Ila 


IC972 


0.00 


0.00 


0.50 


0.50 


0.00 


HaB 


iib/ni 


IC 1297 


0.00 


0.13 


0.87 


0.00 


0.00 


EaB 


nb 


IC 1454 


0.00 


0.21 


0.06 


0.73 


0.00 


HaE 


m 


IC 1747 


0.00 


0.88 


0.02 


0.10 


0.00 


IlaC 


Ila 


IC 2003 


0.00 


0.02 


0.98 


0.00 


0.00 


IlbF 


lib 


IC 2149 


0.00 


0.87 


0.13 


0.00 


0.00 


IlaF 


Ila 


IC2165 


0.00 


0.24 


0.76 


0.00 


0.00 


IlbF 


lib 


IC 2448 


0.00 


0.44 


0.56 


0.00 


0.00 


IlaF 


lib 


IC 2501 


0.02 


0.82 


0.13 


0.03 


0.00 


IlaF 


Ila 


IC 2553 


0.00 


0.88 


0.06 


0.06 


0.00 


IlaC 


Ila 


IC 2621 


0.01 


0.99 


0.00 


0.00 


0.00 


IlaF 


Ila 


IC 3568 


0.00 


0.00 


0.00 


0.67 


0.33 


IVC/IIbB 


III 


1C4191 


0.01 


0.59 


0.40 


0.00 


0.00 


IlbF 


Ila 


IC 4406 


0.00 


0.08 


0.92 


0.00 


0.00 


IlbF 


lib 


IC 4593 


0.00 


0.00 


0.00 


0.09 


0.91 


IVC/IIbB 


IV 


IC 4634 


0.00 


0.00 


0.97 


0.01 


0.02 


IlbF 


lib 


IC 4642 


0.00 


0.00 


0.00 


1.00 


0.00 


(IIa/IIb)C*/(m/IV)B* 


III 


IC 4673 


0.29 


0.71 


0.00 


0.00 


0.00 


VF 


Ila 


IC 4732 


0.00 


0.00 


0.00 


0.00 


1.00 


IIbC/(III/IV)B 


IV 


IC 4776 


0.00 


0.02 


0.97 


0.01 


0.00 


nbF 


nb 


IC 4846 


0.00 


0.00 


0.01 


0.98 


0.01 


nbC/IIIB 


m 


IC 4997 


0.00 


0.00 


0.00 


1.00 


0.00 


(IIa/IIb)C*/(m/IV)B* 


III 


IC5117 


0.01 


0.60 


0.39 


0.00 


0.00 


IlbF 


Ila 


IC 5148-50 


0.00 


0.94 


0.06 


0.00 


0.00 


IB 


Ha 


IC 5217 


0.00 


0.06 


0.00 


0.93 


0.01 


IlaC 


III 


IRAS 15154-525 


1.00 


0.00 


0.00 


0.00 


0.00 


(V/I/IIa)D* 


I 


IRAS 16455-345 


0.03 


0.43 


0.34 


0.16 


0.04 


(V/IIa/IIb/III/IV)D* 


Ila 


J320 


0.00 


0.00 


0.00 


1.00 


0.00 


IIIC/IIbB 


III 


J900 


0.00 


0.01 


0.99 


0.00 


0.00 


IlbF 


lib 


Jnl 


0.00 


1.00 


0.00 


0.00 


0.00 


indef. 


Ila 


Kl-7 


0.00 


0.28 


0.70 


0.02 


0.00 


IlbF 


lib 


K2-1 


0.00 


0.08 


0.92 


0.00 


0.00 


IlbF 


lib 


K3-60 


0.00 


0.99 


0.01 


0.00 


0.00 


IlaD 


Ila 


K3-61 


0.85 


0.15 


0.00 


0.00 


0.00 


IlaD 


I 


K3-66 


0.00 


0.49 


0.03 


0.01 


0.47 


IlaE 


Ila 


K3-67 


0.00 


0.99 


0.01 


0.00 


0.00 


IlaC 


Ila 


K3-68 


0.00 


0.83 


0.13 


0.04 


0.00 


IlaF 


Ila 


K3-70 


1.00 


0.00 


0.00 


0.00 


0.00 


IF 


I 


K5-1 


0.03 


0.81 


0.08 


0.00 


0.08 


(V/IIa)D 


Ila 


K5-3 


0.00 


0.08 


0.57 


0.35 


0.00 


(V/IIb/III/IV)D 


lib 


K5-4 


0.00 


0.15 


0.65 


0.14 


0.06 


(V/IIb/III/IV)D 


lib 


K5-5 


0.00 


0.50 


0.09 


0.41 


0.00 


(V/IIa)D 


Ila 
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Table 6. continued. 



Name 


/j(type I[y) 


/7(type na[y) 


p(type nbly) 


/>(type m\y) 


p(type IVLy) 


Pre-class 


Post-class 


K5-6 


0.00 


0.17 


0.65 


0.16 


0.02 


(V/nb/III/IV)D 


lib 


K5-7 


0.12 


0.88 


0.00 


0.00 


0.00 


(IIa/V)D 


Ila 


K5-9 


0.02 


0.94 


0.01 


0.03 


0.00 


(V/Ila)D 


Ila 


K5-11 


0.17 


0.83 


0.00 


0.00 


0.00 


(V/na)D 


Ila 


K5-12 


0.00 


0.31 


0.51 


0.08 


0.10 


(V/IIb/III/IV)D 


lib 


K5-13 


0.59 


0.41 


0.00 


0.00 


0.00 


(W)D 


I 


K5-14 


0.02 


0.71 


0.23 


0.02 


0.02 


(V/IIa)D 


Ila 


K5-16 


0.19 


0.80 


0.01 


0.00 


0.00 


(na/V)D 


Ila 


K5-17 


0.00 


0.39 


0.35 


0.25 


0.01 


(V/IIa)D 


Ila 


K5-19 


0.72 


0.28 


0.00 


0.00 


0.00 


(IlaA^)D 


I 


K5-20 


0.00 


0.00 


0.09 


0.90 


0.01 


(IlaA^)D 


III 


K6-3 


0.24 


0.75 


0.01 


0.00 


0.00 


(IIa/V)D 


Ila 


K6-5 


1.00 


0.00 


0.00 


0.00 


0.00 


(I/V)D 


I 


K648 


0.00 


0.00 


0.00 


0.00 


1.00 


IVF 


IV 


KFL2 


0.00 


0.02 


0.00 


0.98 


0.00 


VF 


III 


KFL4 


0.00 


0.10 


0.82 


0.07 


0.01 


VP 


lib 


KFL 19 


0.00 


0.00 


0.00 


0.98 


0.02 


IIIC/IIbB 


III 


Lo4 


0.47 


0.53 


0.00 


0.00 


0.00 


ID 


Ila 


Ml-1 


0.00 


0.00 


0.08 


0.02 


0.90 


naB 


IV 


Ml-4 


0.00 


0.55 


0.45 


0.00 


0.00 


IlaF 


Ila 


Ml-5 


0.00 


0.24 


0.76 


0.00 


0.00 


IlbF 


lib 


Ml-6 


0.00 


0.02 


0.97 


0.01 


0.00 


IlbF 


lib 


Ml-7 


0.00 


0.68 


0.32 


0.00 


0.00 


IlaF 


Ila 


Ml-8 


0.93 


0.07 


0.00 


0.00 


0.00 


IF 


I 


Ml-9 


0.00 


0.00 


0.03 


0.97 


0.00 


IlbC/IIIB 


III 


Ml-11 


0.00 


1.00 


0.00 


0.00 


0.00 


IlaF 


Ila 


Ml-12 


0.00 


0.24 


0.28 


0.48 


0.00 


IlbF 


III 


Ml-13 


0.01 


0.97 


0.02 


0.00 


0.00 


IlaF 


Ila 


Ml-14 


0.00 


0.02 


0.98 


0.00 


0.00 


IlbF 


nb 


Ml-16 


0.00 


1.00 


0.00 


0.00 


0.00 


IB 


Ila 


Ml-17 


0.00 


0.24 


0.76 


0.00 


0.00 


IlbF 


lib 


Ml-18 


0.00 


1.00 


0.00 


0.00 


0.00 


IB 


Ila 


Ml-19 


0.00 


0.82 


0.13 


0.05 


0.00 


VF 


Ila 


Ml -20 


0.00 


0.08 


0.78 


0.14 


0.00 


VF 


lib 


Ml -22 


0.00 


1.00 


0.00 


0.00 


0.00 


VF 


Ila 


Ml -25 


0.00 


0.82 


0.18 


0.00 


0.00 


VF 


Ila 


Ml -26 


0.00 


0.24 


0.76 


0.00 


0.00 


VA 


lib 


Ml -27 


0.00 


0.53 


0.47 


0.00 


0.00 


VF 


Ila 


Ml -29 


0.12 


0.87 


0.01 


0.00 


0.00 


VF 


Ila 


Ml-30 


0.00 


0.85 


0.15 


0.00 


0.00 


VF 


Ila 


Ml-31 


1.00 


0.00 


0.00 


0.00 


0.00 


VF 


I 


Ml-34 


0.00 


1.00 


0.00 


0.00 


0.00 


VF 


Ila 


Ml-35 


0.97 


0.03 


0.00 


0.00 


0.00 


VF 


I 


Ml-37 


0.00 


0.02 


0.00 


0.98 


0.00 


VF 


m 


Ml-38 


0.00 


0.00 


0.00 


0.98 


0.02 


VF 


III 


Ml -40 


0.04 


0.96 


0.00 


0.00 


0.00 


VA 


Ila 


Ml -42 


0.00 


1.00 


0.00 


0.00 


0.00 


VF 


Ila 


Ml -44 


0.00 


0.00 


0.00 


0.07 


0.93 


VF 


IV 


Ml -48 


0.00 


1.00 


0.00 


0.00 


0.00 


IB 


Ila 


Ml-50 


0.00 


0.06 


0.80 


0.14 


0.00 


IlbF 


lib 


Ml-51 


0.15 


0.85 


0.00 


0.00 


0.00 


IF 


Ila 


Ml-54 


0.00 


1.00 


0.00 


0.00 


0.00 


IC 


Ila 


Ml-56 


0.00 


0.00 


0.00 


0.17 


0.83 


IlaC 


IV 


M i - J / 


u.u/ 


u.yz 


U.Ui 


f\ f\f\ 


f\ f\f\ 

u.uu 


liar 


Ila 


Ml-58 


0.04 


0.91 


0.05 


0.00 


0.00 


IlaF 


Ila 


Ml-60 


0.00 


0.99 


0.01 


0.00 


0.00 


IlaF 


Ila 


Ml-61 


0.00 


0.54 


0.46 


0.00 


0.00 


IlaF 


Ila 


Ml -67 


0.00 


0.97 


0.01 


0.02 


0.00 


IlaD 


Ila 


Ml-72 


0.00 


0.09 


0.42 


0.49 


0.00 


IlbD 


III 


Ml-74 


0.00 


0.89 


0.11 


0.00 


0.00 


IlaF 


Ila 
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Table 6. continued. 



Name 


/?(type 


pitypc na[y) 


/>(type Ilbly) 


/)(type m\y) 


p(type TV\y) 


Pre-class 


Post-class 


Ml-75 


1.00 


0.00 


0.00 


0.00 


0.00 


IF 


I 


Ml-78 


0.00 


0.86 


0.00 


0.12 


0.02 


IlaC 


Ila 


Ml-79 


0.95 


0.05 


0.00 


0.00 


0.00 


IlaF 


I 


Ml-80 


0.01 


0.97 


0.02 


0.00 


0.00 


IlaF 


Ila 


M2-2 


0.01 


0.77 


0.22 


0.00 


0.00 


IlaF 


Ila 


M2-4 


0.00 


0.62 


0.38 


0.00 


0.00 


IlaF 


Ila 


M2-6 


0.00 


0.01 


0.51 


0.40 


0.08 


(V/IIb/III/IV)D 


lib 


M2-7 


0.00 


0.00 


0.00 


0.01 


0.99 


VF 


IV 


M2-8 


0.00 


0.00 


0.00 


1.00 


0.00 


VF 


III 


M2-9 


0.00 


0.00 


0.01 


0.99 


0.00 


IIIF 


III 


M2-10 


0.00 


0.97 


0.03 


0.00 


0.00 


VF 


Ila 


M2-11 


0.00 


0.00 


0.00 


0.00 


1.00 


VF 


IV 


M2-12 


0.00 


0.00 


0.00 


1.00 


0.00 


VF 


III 


M2-13 


0.00 


0.00 


0.00 


0.00 


1.00 


IlaE 


IV 


M2-15 


0.00 


0.00 


0.00 


0.00 


1.00 


IlaC 


IV 


M2-16 


0.00 


0.99 


0.00 


0.01 


0.00 


VF 


Ila 


M2-18 


0.00 


0.00 


0.23 


0.76 


0.01 


VF 


III 


M2-19 


0.00 


0.15 


0.25 


0.59 


0.01 


VF 


III 


M2-20 


0.00 


0.24 


0.43 


0.33 


0.00 


VF 


lib 


M2-21 


0.00 


0.00 


0.00 


0.72 


0.28 


VF 


III 


M2-22 


0.01 


0.99 


0.00 


0.00 


0.00 


VF 


Ila 


M2-23 


0.00 


0.00 


0.00 


0.00 


1.00 


VF 


IV 


M2-24 


0.00 


0.00 


0.00 


0.00 


1.00 


VF 


IV 


M2-26 


0.00 


0.00 


0.00 


0.38 


0.62 


VF 


IV 


M2-27 


0.00 


0.00 


0.00 


1.00 


0.00 


VF 


III 


M2-29 


0.00 


0.00 


0.00 


0.00 


1.00 


VF 


IV 


M2-30 


0.00 


0.01 


0.08 


0.91 


0.00 


VF 


III 


M2-33 


0.00 


0.00 


0.00 


0.12 


0.88 


VF 


IV 


M2-36 


0.00 


1.00 


0.00 


0.00 


0.00 


VF 


Ila 


M2-38 


0.00 


0.00 


0.00 


0.00 


1.00 


VF 


IV 


M2-39 


0.00 


0.17 


0.83 


0.00 


0.00 


VF 


lib 


M2-42 


0.00 


0.08 


0.00 


0.63 


0.29 


VF 


III 


M2-46 


0.00 


1.00 


0.00 


0.00 


0.00 


IlaF 


Ila 


M2-50 


0.00 


0.02 


0.05 


0.93 


0.00 


IlbC/IIIB 


III 


M2-52 


1.00 


0.00 


0.00 


0.00 


0.00 


IC 


I 


M2-53 


0.00 


0.00 


0.08 


0.88 


0.04 


IlbC/IIIB 


III 


M2-55 


0.02 


0.98 


0.00 


0.00 


0.00 


IlaF 


Ila 


M3-1 


0.00 


0.19 


0.81 


0.00 


0.00 


IlaF 


lib 


M3-2 


0.00 


1.00 


0.00 


0.00 


0.00 


IB 


Ila 


M3-3 


0.00 


1.00 


0.00 


0.00 


0.00 


IC 


Ila 


M3-4 


0.24 


0.76 


0.00 


0.00 


0.00 


IlaF 


Ila 


M3-5 


0.15 


0.85 


0.00 


0.00 


0.00 


IlaF 


Ila 


M3-6 


0.00 


0.00 


0.90 


0.00 


0.10 


IlbF 


lib 


M3-7 


0.00 


0.00 


0.00 


0.02 


0.98 


VF 


IV 


M3-8 


0.00 


0.33 


0.00 


0.67 


0.00 


VF 


III 


M3-9 


0.00 


0.00 


0.04 


0.96 


0.00 


VF 


III 


M3-10 


0.00 


0.07 


0.01 


0.91 


0.01 


VF 


III 


M3-14 


0.38 


0.62 


0.00 


0.00 


0.00 


VF 


Ila 


M3-15 


0.00 


1.00 


0.00 


0.00 


0.00 


VF 


Ila 


M3-16 


0.00 


0.00 


0.12 


0.78 


0.10 


VF 


III 


M3-17 


0.00 


0.11 


0.89 


0.00 


0.00 


VF 


lib 


M3-19 


0.00 


0.00 


0.00 


0.90 


0.10 


VF 


III 


M3-20 


0.00 


0.01 


0.99 


0.00 


0.00 


VF 


lib 


Mj-Zl 


u.uu 


U.JO 


U.OZ 


U.UU 


u.uu 


Vr 


lib 


M3-26 


0.00 


0.00 


0.00 


0.00 


1.00 


VF 


IV 


M3-29 


0.00 


0.95 


0.05 


0.00 


0.00 


IlaF 


Ila 


M3-32 


0.00 


0.00 


0.00 


0.00 


1.00 


VF 


IV 


M3-33 


0.00 


0.00 


0.00 


1.00 


0.00 


IIIF 


III 


M3-35 


0.00 


0.00 


0.00 


0.00 


1.00 


(IIa/IIb)C*/(III/IV)B* 


IV 


M3-36 


0.00 


0.11 


0.71 


0.17 


0.01 


VF 


lib 
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Table 6. continued. 



Name 


/j(type I[y) 


pitypc na[y) 


/>(type nbly) 


/)(type m\y) 


p(type TV\y) 


Pre-class 


Post-class 


M3-37 


0.00 


0.01 


0.05 


0.58 


0.36 


VF 


III 


M3-38 


0.00 


0.53 


0.00 


0.02 


0.45 


VF 


Ila 


M3-42 


0.00 


0.00 


0.00 


0.00 


1.00 


VF 


IV 


M3-43 


0.63 


0.36 


0.01 


0.00 


0.00 


VF 


I 


M3-45 


0.00 


0.12 


0.88 


0.00 


0.00 


VF 


lib 


M3-46 


0.87 


0.13 


0.00 


0.00 


0.00 


VD 


I 


M3-48 


0.00 


0.94 


0.06 


0.00 


0.00 


VD 


Ila 


M3-50 


0.00 


0.00 


0.00 


0.88 


0.12 


VD 


III 


M3-52 


0.00 


1.00 


0.00 


0.00 


0.00 


IE 


Ila 


M3-54 


0.00 


0.00 


1.00 


0.00 


0.00 


IlbF 


lib 


M4-3 


0.00 


0.00 


0.00 


0.00 


1.00 


VF 


IV 


M4-4 


0.00 


0.00 


0.00 


1.00 


0.00 


VF 


III 


M4-6 


0.00 


0.00 


0.00 


0.00 


1.00 


VF 


IV 


M4-7 


0.00 


0.00 


0.00 


0.00 


1.00 


VF 


IV 


M4-8 


0.00 


0.00 


0.00 


1.00 


0.00 


IlbC/IIIB 


III 


M4-14 


0.01 


0.88 


0.11 


0.00 


0.00 


indef. 


Ila 


MaC2-l 


0.00 


0.00 


0.28 


0.48 


0.24 


(IIb/III/IV)D 


III 


Mel-1 


0.00 


1.00 


0.00 


0.00 


0.00 


IlaC 


Ila 


Me2-1 


0.00 


0.00 


0.00 


0.00 


1.00 


rvF 


IV 


Me2-2 


0.00 


1.00 


0.00 


0.00 


0.00 


IE 


Ha 


MyCnlS 


0.01 


0.99 


0.00 


0.00 


0.00 


IlaF 


Ha 


Mzl 


0.02 


0.97 


0.01 


0.00 


0.00 


IlaF 


Ila 


Mz2 


0.01 


0.85 


0.14 


0.00 


0.00 


IlaF 


Ila 


Mz3 


0.03 


0.97 


0.00 


0.00 


0.00 


IC 


Ila 


NGC 40 


0.00 


0.71 


0.29 


0.00 


0.00 


IlbF 


Ila 


NGC 246 


0.00 


0.49 


0.23 


0.27 


0.01 


indef. 


Ila 


NGC 650 


0.20 


0.79 


0.01 


0.00 


0.00 


IlaF 


Ila 


NGC 1360 


0.00 


0.70 


0.29 


0.01 


0.00 


indef. 


Ila 


NGC 1514 


0.00 


0.72 


0.26 


0.02 


0.00 


(na/nb)F* 


Ila 


NGC 1535 


0.00 


0.00 


0.07 


0.85 


0.08 


IIIC/IIbB 


III 


NGC 2022 


0.00 


1.00 


0.00 


0.00 


0.00 


IlaF 


Ila 


NGC 2242 


0.00 


0.00 


0.00 


0.01 


0.99 


IlaB 


IV 


NGC 2346 


0.09 


0.91 


0.00 


0.00 


0.00 


IlaF 


Ila 


NGC 2371-72 


0.10 


0.89 


0.01 


0.00 


0.00 


IlaF 


Ila 


NGC 2392 


0.00 


0.98 


0.02 


0.00 


0.00 


IlaC 


Ila 


NGC 2438 


0.03 


0.94 


0.03 


0.00 


0.00 


IlaF 


Ila 


NGC 2440 


0.20 


0.80 


0.00 


0.00 


0.00 


IlaF 


Ila 


NGC 2452 


0.10 


0.90 


0.00 


0.00 


0.00 


IlaF 


Ila 


NGC 2474-75 


0.00 


1.00 


0.00 


0.00 


0.00 


indef. 


Ila 


NGC 2610 


0.00 


0.94 


0.04 


0.02 


0.00 


IlaC 


Ila 


NGC 2792 


0.01 


0.62 


0.37 


0.00 


0.00 


IlaF 


Ila 


NGC 2818 


0.00 


1.00 


0.00 


0.00 


0.00 


IB 


Ila 


NGC 2867 


0.05 


0.89 


0.06 


0.00 


0.00 


IlaF 


Ila 


NGC 2899 


0.97 


0.03 


0.00 


0.00 


0.00 


IF 


I 


NGC 3132 


0.02 


0.91 


0.07 


0.00 


0.00 


IlaF 


Ila 


NGC 3195 


0.00 


0.74 


0.26 


0.00 


0.00 


IlaF 


Ila 


NGC 3211 


0.00 


0.26 


0.74 


0.00 


0.00 


nbF 


lib 


NGC 3242 


0.00 


0.13 


0.87 


0.00 


0.00 


IlbF 


lib 


NGC 3587 


0.00 


0.27 


0.73 


0.00 


0.00 


IlbF 


lib 


NGC 3699 


0.03 


0.96 


0.01 


0.00 


0.00 


IlaF 


Ila 


NGC 3918 


0.01 


0.81 


0.18 


0.00 


0.00 


IlaF 


Ila 


NGC 4361 


0.00 


0.07 


0.87 


0.06 


0.00 


IlbF 


lib 


NGC 5189 


0.06 


0.92 


0.02 


0.00 


0.00 


IlaF 


Ila 




U.UU 


U.UU 


U.UU 


v.yy 


U.Ui 


11dL./(111/1V )H 


III 
111 


NGC 5315 


0.73 


0.27 


0.00 


0.00 


0.00 


IlaF 


I 


NGC 5873 


0.00 


0.06 


0.82 


0.05 


0.07 


IlaB 


lib 


NGC 5882 


0.00 


0.95 


0.05 


0.00 


0.00 


IlaF 


Ila 


NGC 5979 


0.00 


0.53 


0.00 


0.47 


0.00 


IlaC 


Ila 


NGC 6058 


0.00 


0.00 


0.00 


0.00 


1.00 


IlaB 


IV 


NGC 6153 


0.64 


0.36 


0.00 


0.00 


0.00 


IC 


I 
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Table 6. continued. 



Name 


/j(type I[y) 


pitypc na[y) 


/>(type nbly) 


/)(type m\y) 


p(type TV\y) 


Pre-class 


Post-class 


NGC 6210 


0.00 


0.00 


0.10 


0.90 


0.00 


IlbF 


III 


NGC 6302 


1.00 


0.00 


0.00 


0.00 


0.00 


IF 


I 


NGC 6309 


0.00 


0.14 


0.47 


0.39 


0.00 


IlbF 


lib 


NGC 6326 


0.00 


0.01 


0.20 


0.79 


0.00 


llbC/lIIB 


III 


NGC 6337 


0.00 


0.90 


0.00 


0.10 


0.00 


IlaF 


Ila 


NGC 6369 


0.00 


0.02 


0.00 


0.98 


0.00 


IlaC 


III 


NGC 6439 


0.00 


0.01 


0.00 


0.99 


0.00 


IE 


III 


NGC 6445 


0.98 


0.02 


0.00 


0.00 


0.00 


IF 


I 


NGC 6537 


1.00 


0.00 


0.00 


0.00 


0.00 


IF 


I 


NGC 6543 


0.00 


0.97 


0.01 


0.02 


0.00 


IlaC 


Ila 


NGC 6563 


0.04 


0.95 


0.01 


0.00 


0.00 


VA 


Ila 


NGC 6565 


0.00 


0.98 


0.02 


0.00 


0.00 


VF 


Ila 


NGC 6567 


0.00 


0.00 


0.00 


0.99 


0.01 


IIbC/(III/IV)B 


III 


NGC 6572 


0.01 


0.97 


0.02 


0.00 


0.00 


IlaF 


Ila 


NGC 6578 


0.00 


0.60 


0.40 


0.00 


0.00 


IlbF 


Ila 


NGC 6620 


0.00 


0.98 


0.02 


0.00 


0.00 


VF 


Ila 


NGC 6629 


0.00 


0.00 


1.00 


0.00 


0.00 


VA 


lib 


NGC 6644 


0.00 


0.00 


0.00 


0.00 


1.00 


VF 


IV 


NGC 6720 


0.06 


0.89 


0.05 


0.00 


0.00 


IlaF 


Ila 


NGC 6741 


0.02 


0.94 


0.04 


0.00 


0.00 


IlaF 


Ila 


NGC 6751 


0.01 


0.98 


0.01 


0.00 


0.00 


IlaC 


Ila 


NGC 6765 


0.00 


0.15 


0.00 


0.84 


0.01 


IlaC 


III 


NGC 6778 


0.97 


0.03 


0.00 


0.00 


0.00 


IC 


I 


NGC 6781 


0.40 


0.60 


0.00 


0.00 


0.00 


IF 


Ila 


NGC 6790 


0.00 


0.32 


0.68 


0.00 


0.00 


IlbF 


lib 


NGC 6803 


0.19 


0.81 


0.00 


0.00 


0.00 


IF 


Ila 


NGC 6804 


0.64 


0.35 


0.01 


0.00 


0.00 


(I/IIa)F* 


I 


NGC 6807 


0.00 


0.00 


0.00 


0.00 


1.00 


IIbC/(III/IV)B 


IV 


NGC 6818 


0.00 


0.47 


0.45 


0.08 


0.00 


IlaF 


Ila 


NGC 6826 


0.00 


0.03 


0.97 


0.00 


0.00 


IlbF 


lib 


NGC 6833 


0.00 


0.02 


0.00 


0.00 


0.98 


IlaC 


IV 


NGC 6853 


0.09 


0.91 


0.00 


0.00 


0.00 


IlaC 


Ila 


NGC 6879 


0.00 


0.23 


0.72 


0.05 


0.00 


IlbF 


lib 


NGC 6881 


0.03 


0.96 


0.01 


0.00 


0.00 


IlaF 


Ila 


NGC 6884 


0.00 


0.94 


0.06 


0.00 


0.00 


IlaF 


Ila 


NGC 6886 


0.00 


1.00 


0.00 


0.00 


0.00 


IlaC 


Ila 


NGC 6891 


0.00 


0.03 


0.97 


0.00 


0.00 


IlbF 


lib 


NGC 6894 


0.00 


0.99 


0.00 


0.00 


0.01 


IlaC 


Ila 


NGC 6905 


0.01 


0.34 


0.65 


0.00 


0.00 


nbF 


lib 


NGC 7008 


0.09 


0.91 


0.00 


0.00 


0.00 


IC 


Ila 


NGC 7009 


0.00 


0.96 


0.03 


0.01 


0.00 


IlaC 


lla 


NGC 7026 


0.02 


0.98 


0.00 


0.00 


0.00 


naF 


Ila 


NGC 7027 


0.00 


0.91 


0.09 


0.00 


0.00 


IlaF 


Ila 


NGC 7293 


0.79 


0.21 


0.00 


0.00 


0.00 


IF 


I 


NGC 7354 


0.01 


0.94 


0.05 


0.00 


0.00 


IlaF 


Ila 


NGC 7662 


0.00 


0.07 


0.93 


0.00 


0.00 


IlbF 


lib 


PBl 


0.00 


0.00 


0.49 


0.48 


0.03 


(IIb/III/IV)D 


lib 


PB3 


0.05 


0.90 


0.05 


0.00 


0.00 


IlaD 


Ila 


PB4 


0.06 


0.76 


0.13 


0.05 


0.00 


(l/lIa)D* 


Ila 


PB6 


0.00 


1.00 


0.00 


0.00 


0.00 


IE 


Ila 


PB8 


0.00 


0.80 


0.14 


0.06 


0.00 


(IIa/IIb)F* 


Ila 


PC12 


0.00 


0.00 


0.00 


0.82 


0.18 


IVC/IIbB 


III 


PC14 


0.00 


0.01 


0.03 


0.96 


0.00 


IlbC/IIIB 


III 


Dal 1 

rei-i 


u.uu 


U.jU 


U.4/ 


U.U3 


u.uu 


HUE? 

llbr 


lla 


Pel-12 


0.00 


0.02 


0.15 


0.83 


0.00 


VF 


III 


Pel-13 


0.00 


0.00 


0.00 


0.00 


1.00 


IIaE*/IIIF*/IVB* 


IV 


Pel-17 


0.17 


0.83 


0.00 


0.00 


0.00 


IlaC 


Ila 


Pel-18 


0.30 


0.70 


0.00 


0.00 


0.00 


ID 


Ila 


Pe2-7 


0.00 


0.28 


0.58 


0.11 


0.03 


(IIb/m/IV)D 


lib 


Pe2-ll 


0.00 


0.00 


0.00 


1.00 


0.00 


VD 


III 
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Table 6. continued. 



Name 


/j(type I[y) 


/7(type na[y) 


/>(type nbly) 


/)(type m\y) 


p(type IVLy) 


Pre-class 


Post-class 


Pe2-14 


0.00 


0.88 


0.00 


0.12 


0.00 


IlaF 


Ila 


PNG006.0-41.9 


0.00 


0.00 


0.00 


0.00 


1.00 


IlaE 


IV 


PNG135.9+55.9 


0.00 


0.29 


0.31 


0.34 


0.06 


(na/nb/ni/iv)D* 


III 


PNG243.8-37.1 


0.00 


0.00 


0.00 


0.00 


1.00 


IlaE 


IV 


PWl 


0.95 


0.05 


0.00 


0.00 


0.00 


indef. 


I 


SA2-21 


0.00 


1.00 


0.00 


0.00 


0.00 


IE 


Ila 


SBOl 


0.00 


0.00 


0.37 


0.60 


0.03 


VE/(IIb/III/IV)D 


III 


SB02 


0.64 


0.36 


0.00 


0.00 


0.00 


VE/IIaD 


I 


SB03 


1.00 


0.00 


0.00 


0.00 


0.00 


VE/ID 


I 


SB04 


0.22 


0.72 


0.00 


0.00 


0.06 


VE/IIaD 


Ila 


SB06 


0.91 


0.09 


0.00 


0.00 


0.00 


VE/ID 


I 


SB12 


0.05 


0.47 


0.31 


0.15 


0.02 


(V/IIa/IIb/III/IV)D* 


Ila 


SB15 


0.00 


0.00 


0.00 


0.00 


1.00 


(V/IIb/III/IV)D 


IV 


SB17 


0.25 


0.74 


0.01 


0.00 


0.00 


IlaD 


Ila 


SB18 


0.00 


0.00 


0.02 


0.04 


0.94 


(IIb/in/IV)D 


IV 


SB19 


0.00 


0.05 


0.01 


0.94 


0.00 


IlaD 


III 


SB20 


0.00 


0.00 


0.04 


0.96 


0.00 


IlaD 


III 


SB21 


0.00 


0.00 


0.19 


0.79 


0.02 


(IIb/III/IV)D 


III 


SB24 


0.00 


0.00 


0.11 


0.23 


0.66 


(nb/III/IV)D 


IV 


SB25 


0.94 


0.06 


0.00 


0.00 


0.00 


ID 


I 


SB26 


0.00 


0.07 


0.61 


0.19 


0.13 


(IIb/III/IV)D 


nb 


SB28 


0.04 


0.46 


0.32 


0.15 


0.03 


(IIa/IIb/III/IV)D* 


Ila 


SB30 


0.00 


0.04 


0.62 


0.23 


0.11 


(IIb/III/IV)D 


lib 


SB31 


0.00 


0.05 


0.66 


0.25 


0.04 


(nb/III/IV)D 


lib 


SB32 


0.11 


0.86 


0.03 


0.00 


0.00 


IlaD 


Ila 


SB33 


0.01 


0.90 


0.00 


0.00 


0.09 


VE/IIaD 


Ila 


SB34 


0.56 


0.39 


0.02 


0.03 


0.00 


VE*/(I/IIa)D* 


I 


SB35 


0.00 


0.00 


0.37 


0.45 


0.18 


(V/IIb/III/IV)D 


III 


SB37 


0.73 


0.24 


0.01 


0.02 


0.00 


(I/IIa/V)D* 


I 


SB38 


0.07 


0.86 


0.06 


0.01 


0.00 


(IIa/V)D 


Ila 


SB42 


0.00 


0.00 


0.08 


0.92 


0.00 


VE/(IIb/III/IV)D 


III 


SB44 


0.00 


0.01 


0.13 


0.86 


0.00 


VE/IIaD 


III 


SB50 


0.03 


0.85 


0.08 


0.04 


0.00 


VE/IIaD 


Ila 


SB52 


0.00 


0.00 


0.11 


0.89 


0.00 


(IIa/V)D 


III 


SB53 


0.45 


0.55 


0.00 


0.00 


0.00 


(I/V)D 


Ila 


SB55 


0.00 


0.01 


0.02 


0.97 


0.00 


(IIa/V)D 


III 


Sh2-71 


1.00 


0.00 


0.00 


0.00 


0.00 


indef. 


I 


SkW13-2 


0.00 


0.00 


0.11 


0.89 


0.00 


(V/IIb/III/IV)D 


III 


Snl 


0.00 


0.00 


0.00 


0.00 


1.00 


rvF 


IV 


SwStl 


0.00 


0.01 


0.99 


0.00 


0.00 


VA 


lib 


Tel 


0.00 


0.01 


0.98 


0.01 


0.00 


IlbF 


lib 


Th2-A 


0.00 


0.94 


0.06 


0.00 


0.00 


IlaF 


Ila 


Th3-13 


0.00 


1.00 


0.00 


0.00 


0.00 


VF 


Ila 


Th3-15 


0.00 


0.85 


0.14 


0.01 


0.00 


VD 


Ila 


Th3-34 


1.00 


0.00 


0.00 


0.00 


0.00 


(I/V)D 


I 


Th4-2 


0.00 


1.00 


0.00 


0.00 


0.00 


VF 


Ila 


Th4-7 


0.00 


0.15 


0.05 


0.80 


0.00 


VF 


III 


Th4-10 


0.03 


0.79 


0.15 


0.01 


0.02 


IlaD 


Ila 


Vdl-1 


0.00 


0.00 


0.39 


0.57 


0.04 


(IIb/III/IV)D 


III 


Vyi-l 


u.uu 


A 

0.22 


0. lb 


A AO 
0.02 


A AA 

0.00 


1B*/Ilar* 


lib 


Vyl-2 


0.00 


0.00 


0.00 


0.00 


1.00 


IVF 


IV 


Vy2-1 


0.00 


1.00 


0.00 


0.00 


0.00 


VF 


Ila 


Vy2-2 


0.00 


0.00 


0.00 


0.99 


0.01 


IIbC/(III/IV)B 


III 


Y-C2-5 


0.00 


0.02 


0.07 


0.91 


0.00 


(IIa/IIb)D* 


III 


Y-C2-32 


0.00 


0.00 


0.44 


0.43 


0.13 


(IIb/III/IV)D 


lib 


Ym29 


0.90 


0.10 


0.00 


0.00 


0.00 


indef. 


I 



